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1 Introduction 



Les accelerateurs terrestres peuvent fournir des informations jusqu'a des echelles d'en- 
viron 1 TeV (10 12 eV). L'echelle de Grande Unification (GUT, le point de convergence 
des couplages electrofaible et fort) est aux environs de 10 15 GeV=10 24 eV. L'echelle de 
Planck, a 10 19 GeV=10 28 eV, correspond a la transition vers une gravitation quantique. 
On pourra difficilement obtenir des informations dans ces domaines au moyen des acce- 
lerateurs terrestres. La Cosmologie nous offre un excellent laboratoire pour les etudes a 
haute energie. 

La physique des astroparticules est le trait d'union entre l'Astrophysique, la Cosmo- 
logie et la physique de particules. 

La liste de problemes qui se posent au physicien des particules et a l'astrophysicien 
est tres grande : 

• Le probleme de la matiere obscure et du role de la constante cosmologique ("energie 
noire") . 

• La difficulte que Ton rencontre a concilier l'homogeneite de la radiation de fond a 
2.725 K et la distribution desordonnee de la matiere visible. 

• L'existence de rayons cosmiques de tres haute energie (> 10 20 eV). 

• L'absence, jusqu'a aujourd'hui, de tout signal clair en faveur d'une physique au-dela 
du Modele Standard, que Ton sait incomplet. 

etc... 

La solution de ces problemes, qui demande des efforts materiels et intellectuels consi- 
derables, est fort probablement indispensable a la comprehension de la realite et nous 
amenera certainement a la decouverte de nouveaux mondes (on peut toujours rever!). 
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1.1 Astroparticules ; de la production a l'observation 

Les astroparticules (AP) sont les signaux qui nous sont envoyes des astres et du Cos- 
mos : photons, hadrons, leptons, ondes gravitationnelles,... 

Les astroparticules sont produites dans une source extra terrestre et se propagent vers 
la Terre ou elles sont alors observees par le physicien (figure 1.1). 



Source: 
production, acceleration, 
self- absorption, ... 




Fig. 1.1 - Production, propagation et detection des astroparticules. 
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La production peut etre un due a un phenomene d'acceleration ou a la disintegration 
d'un objet plus ou moms exotique. 

La source peut presenter un phenomene d'autoabsorption. 

Pendant son voyage vers la terre, TAP peut se desintegrer, interagir avec le milieu 
intergalactique, interstellaire, interplanetaire et l'atmosphere. Les champs magnetiques 
modifient la trajectoire des particules chargees. L'interaction avec la radiation e.m. (la 
lumiere des etoiles, la "Cosmic Microwave Background Radiation (CMBR)", les champs 
magnetiques,...) degrade l'information initiate en redistribuant l'energie. L'AP etant arri- 
vee sur terre, 1' observation se fait par : 

• Des detecteurs installes sur navettes spatiales, ayant des volumes typiques de l'ordre 
du m 3 . 

• Des detecteurs au sol ; par exemple des reseaux de scintillateurs couvrant des surfaces 
de l'ordre de 10-100 km 2 . Un autre exemple est donne par les reseaux (Flye's Eye) qui 
utilisent l'atmosphere comme milieu de detection (fluorescence ou lumiere Cherenkov). 

• Des detecteurs de grand volume, dans le sous-sol, pour la detection des neutrinos. 

Si Ton s'interesse a l'extreme du spectre hadronique (E > 10 20 eV), le taux est de 
l'ordre d'un evenement par km 2 , par steradian et par siecle. Seul un detecteur de tres 
grande surface permet l'observation de ces evenements rares. 

Q. : calculer les parametres geometriques d'un detecteur qui utilise l'atmosphere 
comme milieu de conversion et qui utilise la lumiere de fluorescence ou Cherenkov pour 
la detection. On cherche a obtenir au moins 1 evenement /annee. 

Au niveau de la mer, apres 10000 kg/m 2 d'air, on observe au sol essentiellement la 
radiation e.m. non absorbee (onde radio, lumiere), les muons des gerbes hadroniques, 
et les neutrinos. La reconstitution de la composition d'origine (avant interaction avec 
l'atmosphere) n'est pas toujours facile. L'utilisation d'instruments dans des satellites est 
done tres importante, en tout cas pour l'identification des evenements pas trop rares. 

La composition d'origine est donne dans la table, qui oublie les neutrinos, sans compter 
d'autres particules exotiques, pas encore observees. 

98% hadrons hadrons (H/He/autres=87/12/l) 

2% electrons 

0.01% gamma (E > 1 MeV) 

Tab. 1.1 Composition d'origine de la radiation cosmique 
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1.2 Quelques formules et parametres 

1 UA (unite astronomique) 
1 pc (parsec .parallaxe second) 
1 ly (annee lumiere) 
1 masse solaire 
1 rayon solaire 
Luminosite du soleil 
Vitesse du soleil dans la galaxie 

Tab. 1.2 Systeme solaire 



= 1.496 10 11 m 
= 3.0856 10 16 m 
= 9.4605 10 15 m 
= 1.989 1 30 kg 
= 6.9598 10 8 m 
= 3.90 10 26 W 
= 250 km/s 



Donnees sur notre galaxie 



Epaisseur au centre 
Diametre du disque 
Halo 

Distance du soleil au centre de la galaxie 
Epaisseur de la galaxie a la position du soleil 
Hauteur du soleil au-dessus du plan galactique 
Pole Galactique 
Centre de la galaxie 
Masse visible 
Masse (avec halo HI) 
Densite 



= 5 kpc 
= 30 kpc 

= 3x10 kpc (amas globulaires,...) 
= 10 kpc 
= lkpc 
= 10 pc 

= (a = 12h49m, 5 = +27° 24') 
= (a = 17h42m, 5 = -28° 57') 
= 2 10 11 masses solaires 
= 6 10 11 masses solaires 
= environ 11~ 24 g cm -3 



Tab. 1.3 Notre galaxie 



T fond cosmique 
Champ B intergalactique 

galactique 

interplanetaire 

terrestre 

Energie liberee par une Super Nova (SN) 



2.725 ± 0.001 K 
0(1 nG) 
0(1 /iG) 
w 10 fiG 

0.45 G (sous notre latitude) 
10 46 J 



Tab. 1.4 - Varia 



Rayon de courbure pour une particule d'energie totale E, de quantite de mouvement 
p, de masse m et de charge q=Ze, dans un champ B : 



r = p/(qB) (1.1) 



^-= 200^1^ (1.2) 
1UA ITeV B y ' 
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ou 

1 E 

r[m] = 3 x 10 9 7— - ou 7 = — - (1.3) 

1 J 'B/10- 9 Tesla mc 2 K ' 

Rayon de Schwarzschild = 2GM/c 2 (1.4) 

, , aft) 

Parametre de Hubble H(t) = -+4 (1-5) 

v 7 a(f) v 7 

if (maintenant) = H = h x km s~ 1 Mpc~ 1 (1.6) 

/i = 71 ±0.07 ±(+1.15, -0.95) (1.7) 

"red shift" : (1 + z) = A o6s /A (1.8) 

1.3 Les unites naturelles et les sections efficaces 

Si Ton pose c = h = k B = 1 alors 

longueur =temps= 1 /masse 

lGeV" 2 = 0.389 10- 27 cm 2 (1.9) 

Rayon de Schwarzschild = 2GM (1.10) 

« em (0GeV) = £- « ^ (1.11) 

a em (100GeV) « ^ (1.12) 

charge faible = p^y = e/ sin 6V sin 2 0.23 (1-13) 

^wa em x 0.23 (1.14) 

4tt v ' 

a strong (10GeV) « 0.3 a strong (100GeV) « 0.11 (1.15) 
L'energie dans le cm. est ^/s. 



Sections efficaces a+h -> o + d 




M 






Le propagateur pour une particule de masse M donne une contribution S _ 1 M2 a l'am- 
plitude. 

La section efficace a la forme : 



a = const x p x 



a 



(s-M 2 y 



(1.16) 



a est la constante de l'intensite de l'interaction en jeu ; 

p est le facteur d'espace de phase qui mesure l'ensemble des configurations disponibles 
dans l'etat final. Pour une energie plus grande que la valeur des masses des particules 
dans l'etat final, p croit comme s. Done 



sa 



(J r*j 



(s-M 2 ) 2 



(1.17) 



Si M 2 < s alors a ~ — 



(1.18) 



Si M > s alors a 



sa 
M 1 



(1.19) 



Si Ton se place a des valeurs de s proches de M 2 il faut modifier la section efficace 
ci-dessus en introduisant une forme de Breit-Wigner. Le propagateur devient s _ M 2 +iTM 
avec T la largeur de la resonance. 

Q. : Pour le processus V e e~ — > V^pT calculer Tenergie du neutrino pour que Ton tombe 
sur la masse du W. L' electron est au repos (Reponse : 6.3 PeV). 

Exemples : 



a(e + e — > p + p ) = 



Anal 
3s 



a grand s 



(1.20) 



a(v e e -> v^p ) = 



(1.21) 



a( 11 ^e+e-) = -^(l-v 



2m 2 e 



+ 2v(v z - 2) 



(1.22) 
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avec v la vitesse d'un des leptons dans le cm. : v = \J\ — Ami/ s 



Q. : Montrer que pres du seuil <r(77 — > e + e ) = 
tandis que pour s ^> (2m e ) 2 (7(77 — > e + e~) = 



hi 



Annihilation 



(Tie e 



77) 



2vm1 



(W) 



2w 



1 - v 



(1.23) 



10" l4 b- 



!^ ! 1 1 I 



r ■ j r 




Fig. 1.2 - Sections efiicaces pour rannihilation, la diffusion Compton et la creation de 
paire. 



Diffusion Compton : 



a {^e — > 7e) 



7TCC 



v 3 ml 



[l-v) 



+ { v * + 2v-2 In - 

1 + 1> \1 — i; y 



2w 3 (l + 2u) 
(1 + ^) 2 



(1.24) 



avec v la vitesse de l'electron dans le cm. : v — {s — m 2 e )/{s + m 2 e ) 

Q. : montrer que pour des petites energies, on trouve la formule de Thomson : 

<r(7e — * 7e) = 8 "'" 



*2 



3m2 
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La figure 1.2 montre ces sections efficaces en fonction de la variable y qui donne 
l'energie disponible dans le cm. : 




^/Smin = 2m e pour la creation de paire et pour l'annihilation, y / s^ m " = m e pour la 
diffusion Compton. 

Interaction forte : 

Dans ce cas, la structure interne des hadrons joue un role. 

Dans une diffusion nucleon-nucleon a basse energie, on peut considerer que l'echange 
d'impulsion se fait par un pion et 

a(pp -> X) ~ -At « 5 = 0.4 x 10" 23 cm 2 (1.25) 

Kyy ! ml 0.1 2 GeV 2 V ; 

A partir du moment ou l'interaction est capable de resoudre la structure interne de l'ha- 
dron, on doit faire intervenir les fonctions de structure de ce dernier. Les fonctions de 
structure indiquent la probability pour un parton de transporter une certaine fraction de 
la quantite de mouvement totale de l'hadron. 

Certaines section efficaces sont tabulees (Particle Data Group, PDG). 



1.4 Gerbes hadroniques et e.m. dans l'atmosphere 

Lors du developpement d'une gerbe e.m. (E>0.1 GeV), la distance a laquelle on trouve 
le nombre le plus eleve de particules est 

E , 
In 2 e c 

ou X = longueur de radiation et e c = energie critique 
le nombre de particules au maximum de la gerbe vaut 

Xmax = - (1-27) 

Dans le cas des gerbes hadroniques on a aussi les proportionnalites 

ln£, N max (xE (1.28) 

Une regie valable pour les tres hautes energies est que l'energie de la particule initiale 
peut etre obtenue (a ±25%) en comptant 1.4 GeV par particule au maximum de la gerbe. 



1.5 Les distributions a grande echelle 

Le decouplage entre la radiation et la matiere a eu lieu quand Tunivers avait une taille 
d'environ 1/1000 de celle d'aujourd'hui (z = 1100). L'age etait d'environ td = 300000 
annees. 
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La radiation de fond CMBR a T=2.728± 0.002 K represente revolution du corps noir 
de td a aujourd'hui. Le spectre est celui d'un corps noir avec une densite d'energie 

£rad = crT 4 = 4.2 10" 14 J m- 3 = 2.62 10 5 eV m" 3 (1.29) 

et une densite d'environ 420 photons/cm 3 . 

Le Cosmic Background Explorer (COBE) a effectue des mesures a 31.5, 53 et 90 GHz. 
La mesure brute montre que le fond constitue une espece de repere absolu dans lequel la 
terre se deplace a une vitesse de 200-250 km/s. Une fois corrige l'effet du mouvement de 
la terre et les differentes sources de bruit etant soustraites, le fond restant apparait tres 
uniforme. Si Ton considere des regions angulaires de 10°, les fluctuations sont de 29±l/x K. 
Converti en termes d'intensite, cela donne 

AI/I = 10~ 5 (1.30) 

Cela est montre a la figure 1.3 (voir site web de COBE). Plus recemment le satellite 
WMAP a effectue des observations avec une resolution plus poussee (http :/ /map. gsfc.nasa.gov/) 
a 23, 33, 41, 61 et 94 GHz (figure 1.4). 

Done la radiation presente une homogeneite remarquable. On ne peut pas en dire 
autant pour la matiere "visible". En utilisant la loi de Hubble, on peut estimer les distances 
des galaxies. La figure 1.5 montre la distribution d'environ 100000 galaxies de l'Univers 
proche. II s'agit de galaxies avec des vitesses d'eloignement jusqu'a 15000 km/s. Nous 
sommes au centre de la figure dans l'Amas Local qui, a son tour, se trouve dans le 
Superamas Local. Ce dernier est un ensemble de galaxies centrees sur l'amas Virgo qui se 
trouve a 15-20 Mpc de chez nous. La figure montre aussi la "Grande Muraille". 

La figure 1.6 montre une image de l'Univers local dans laquelle on ne considere que 
les structures les plus importantes. 

Un des problemes de la Cosmologie est de concilier l'homogeneite de la distribution 
du CMBR avec l'irregularite de la matiere visible. Si Ton se limite au domaine des sources 
radio, on selectionne quasi exclusivement les radiogalaxies et les quasars, tous des objets 
tres eloignes. Dans ce cas la distribution est tres uniforme, comme le montre la figure 1.7. 
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Fig. 1.3 - Distribution de l'intensite de la CMBR. En haut : resultat brut de la mesure 
par le "Differential Microwave Radiometer" de COBE. Au milieu : resultat corrige par 
soustraction de l'effet du mouvement de la terre et du bruit de fond du aux contributions 
thermiques (poussieres) et aux gaz chauds (Compton). La radiation synchrotron de notre 
galaxie se manifeste par une bande uniforme traversant la figure de part en part. En bas : 
resultat final, apres soustraction de la contribution de cette radiation synchroton. 
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Fig. 1.4 - Distribution de l'intensite de la CMBR par WMAP. 



Structure autour du Superamas Local 

t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




Perse us -Pisces 



I I I I I I I I I J I I I I I I I I I 

-100 100 

X (Mpc) 

Fig. 1.5 - Vue en coupe de l'univers local a partir d'un echantillon de plus de 100 000 
galaxies. 
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Fig. 1.6 - Structure locale. Nous sommes au centre, sur la petite croix, Taxe X est dirige 
vers le centre de notre Galaxie et l'axe Y est parallele a la vitesse de rotation de celle-ci. 
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Fig. 1.7 - Distribution de 31000 sources radio. Les trous correspondent a des sources 
intenses. Mis a part un leger exces en correspondance du plan galactique, la distribution 
de ces objets tres eloignes est uniforme. 
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2 Les rayons cosmiques 



Les rayons cosmiques (RC) ont ete decouverts en 1912 par Victor Hess par l'etude de 
leur effet ionisant. L'etude se faisait avec des detecteurs emportes dans des ballons. 

Hess a partage le prix Nobel avec C. D. A. Anderson qui avait decouvert le positon 
dans les rayons cosmiques observes dans des emulsions photographiques. 

A partir de 1 TeV les particules du rayonnement cosmique arrivent sur la terre de 
fagon pratiquement isotrope. A plus faible energie, les champs magnetiques du systeme 
solaire affectent la trajectoire de ces particules. 

Q. : montrer que pour p< 10 19 eV les rayons cosmiques charges (p. ex. des protons, 
Z=l) sont confines dans notre galaxie. 

Par "rayons cosmiques" on entend souvent la composante hadronique du flux de par- 
ticules qui entrent dans l'atmosphere terrestre. Dans ce chapitre on va aussi mentionner 
la composante electronique. 

2.1 Spectre d'energie et composition de la composante hadro- 
nique a basse energie 

Le spectre jusqu'a 10 6 MeV/nucleon est represents pour H, He, C, Fe a la figure 2.1. 
On voit que le gros de la statistique se situe entre 10 1 et 10 4 MeV. 




10 1 10= 10 3 10* 10 s 10 s 10 T 

Kinetic onarav/MsV nuclaon-' 



Fig. 2.1 - Spectres mesures en dehors de l'atmosphere pour les noyaux de H, He, C et 
Fe. 

Tous les elements chimiques ont ete observes dans les RC avec des abondances relatives 
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similaires a celles du systeme solaire et de la region interstellaire locale. La figure 2.2 
compare ces deux distributions, en fonction de l'element chimique. Les deux distributions 
sont arbitrairement normalisees au Si (100). 

On observe : 

1) un taux inferieur de protons et d' Helium dans le rayonnement cosmique a celui du 
systeme solaire 

2) un taux superieur pour le groupe Li, Be, B et pour le groupe Sc, Ti, V, Cr, Mn. 

L'observation 1) peut etre le reflet d'une composition a la source plus riche en noyaux 
lourds, mais aussi le fait que ces derniers sont plus facilement ionisables, done plus fa- 
cilement prets a etre acceleres. L'observation 2) est probablement due au phenomene de 
spallation. A partir de C et O, la spallation enrichit le premier groupe, tandis que le Fe 
enrichit le second. 




Fig. 2.2 - Composition chimique du systeme solaire (en continu) et du rayonnement 
cosmique (en pointille). 
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2.2 Les energies ultra- hautes (UHE) 

La figure 2.3 montre le spectre d'energie differentiel a partir de 0.1 GeV. 

T 10< 



> 



^io z - 



>§ -r 

-10 



10 



■4_ 



10" 



10 



10 



10 



10 



10 



10 



10 



■10 
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16 _ 



19- 



22- 



25- 
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Fluxes of Cosmic Rays 



(1 particle per m 2 -second) 



-A 



Knee 

(1 particle per m 2 -year) 



Ankle 

(1 particle per km^-year) 




' 



J i i i 



mwl 



ul I 



' i 



10 9 10 10 10 11 10 12 10 13 10 14 10 15 10 16 10 17 10 18 10 19 10 20 10 21 

Energy (eV) 



Fig. 2.3 - Compilation de mesures du spectre d'energie differentiel des rayons cosmiques. 
La ligne pointillee correspond a a = 3 (voir eq. 2.1 et 2.2). Les points de changement 
d'exposant sont indiques, ainsi que les flux approximatifs par steradian. 
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A partir de 10 GeV/nucleon le spectre plonge avec une loi de puissance 



(2 ,) 



La valeur du parametre est : 




E < 10 16 eV 

a = { 3.0 10 16 < E < 10 18 eV (2.2) 
E > 10 18 eV 

Les points de changement d'indice sont appeles le genou ("knee") et la cheville ("ankle"). 

Une compilation plus recente est montree a la figure 2.4, ou Ton a aplati la distribution 
en multipliant l'ordonnee par E~ 2 7 . A la rigueur, une quatrieme region d'indice a = 2.5 
est visible pres du genou. 

En octobre 1991, le detecteur Fly's Eye a enregistre un evenement qui a ete estime a 
3.01^4 x 10 20 eV. D'autres evenements de ce type ont ete observes par les observatoires 
d'Akeno (Japon) et d'Yakutsk (Russie). 

Peut-on determiner la source des RC ? Si Ton se concentre sur ces energies extremes, 
on peut supposer que l'effet des champs magnetiques est petit. Si Ton considere l'effet des 
champs intergalactiques, la deviation angulaire depend de la distance et du champ. Dans 
les unites appropriees au probleme : 

a -A B 

_ iMpc^G (2 3) 



0.1 



3xl0 20 eV 



Les speculations sur la valeur de B vont de 10 7 a 10 12 Gauss. Ce qui fait que la precision 
sur une centaine de Mpc est comprise entre une fraction de degre et l'infini (ou presque). 

Q. : chercher dans la litterature les methodes utilisees pour la mesure de B. 

D'autre part, si Ton considere des sources ponctuelles a l'interieur de la galaxie, les 
particules avec E > 10 19 eV ont des rayons de courbure plus grands que la taille de la 
galaxie et leur trajectoire devrait pointer vers l'origine des particules. Toutefois aucune 
source n'a pu etre identifiee jusqu'a present. En particulier, aucune correlation avec le plan 
galactique n'a ete mise en evidence, ce qui suggere fortement une origine extragalactique 
pour les rayons cosmiques UHE. 

On peut aussi observer que l'anisotropie dans la distribution de la direction d'arrivee 
des cosmiques est tres petite. Dans la fenetre 10 12 — 10 14 eV l'analyse par ajustment de 
fonctions harmoniques donne des deviations plus petites que 10~ 3 . 



2.2.1 Detection des rayons cosmiques UHE 

L'observation de rayons cosmiques UHE (E >PeV) demande des moyens considerables 
car le taux se situe au niveau d'un evenement par km 2 par siecle, des que Ton depasse 10 20 
eV. La mesure de l'energie necessite une forme de calorimetrie. De plus, il est souhaitable 
de pouvoir determiner la direction de la particule primaire. 
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Fig. 2.4 - Spectre d'energie differentiel des rayons cosmiques. Pour aplatir la distribution, 
on a multiplie le taux par E~ 2 - 7 (en haut) et E~ 3 (le zoom en bas). 



En pratique, on peut utiliser deux methodes : 

1) on compte les particules residuelles de la fin de la gerbe, c. a d. les muons qui 
sont produits par la disintegration des pions. On utilise par exemple une matrice de 
scintillateurs, comme le reseau AGASA au Japon qui couvre 100 km 2 ; 
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2) on utilise la lumiere de fluorescence de l'azote de Pair excite par les particules de 
la gerbe. C'est le cas du detecteur Fly's Eye, dans l'Utah (USA). On peut aussi utiliser 
la lumiere Cherenkov produite par ces particules. Jusqu'a present cette technique n'a ete 
appliquee qu'a la detection de gamma cosmiques. 

Dans les deux cas, l'utilisation d'un programme Monte Carlo pour la modelisation 
de la gerbe est indispensable. De ces simulations on verifie la regie vue au paragraphe 
1.4 : l'energie de la particule initiale peut etre obtenue (a ±25%) en comptant 1.4 GeV 
par particule, au maximum de la gerbe. Done le probleme se ramene (pour ainsi dire) a 
compter le nombre de particules au maximum de la gerbe. Une gerbe de 10 20 eV comporte 
de l'ordre de 10 11 particules, ce qui rend la simulation assez difficile ! 

Le temps relatif d'arrivee des particules de la gerbe permet d'estimer la direction 
d'arrivee de la particule originelle (figure 2.5). D'autres moyens seront discutes plus loin 
(IACT, Imaging Air Cherenkov Telescopes). 




reseau de detecteurs 



Fig. 2.5 - Le front de la gerbe est presque plan et tres compact (dispersion inferieure 
ou de l'ordre d'une ns). Les differences de temps de detection peuvent etre utilisees pour 
estimer l'angle d'entree de la particule primaire dans l'atmosphere. 
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Fig. 2.6 - Anisotropic observee par le detecteur AGASA. 



Le Akeno Giant Air Shower Array (AGASA) est situe dans un vignoble au Ja- 
pon. II couvre 100 km 2 , avec plus de 111 scintillateurs a la surface et 27 en profon- 
deur (pour les muons). Pour plus de details, voir le site web http ://www-akeno.icrr.u- 
tokyo.ac.jp/AGASA/. 

La figure 2.6 montre l'anisotropie observee par AGASA. Les grands carres fonces sont 
les sources avec plus de 10 20 eV et les points clairs les sources entre 4 10 19 et 10 20 eV. 

Le detecteur HEGRA 

Le detecteur "High Energy Gamma Ray Astronomy" (voir figure 2.7) est situe aux 
Canaries, il est plus precisement consacre aux gamma. II comprend actuellement : 

105 stations de scintillateurs pour les particules e.m., avec un seuil de 20 TeV 
49 stations pour la detection de la lumiere Cherenkov, avec un seuil 15 TeV 
17 stations pour les muons et la calorimetrie (telescopes Cherenkov). 

Voir le site web http ://wpos6. physik.uni-wuppertal.de :8080/ 
Le detecteur HiRes Fly's Eye 

Ce detecteur est une amelioration de l'ancien Fly's Eye (voir figure 2.8). La figure 2.9 
montre les tonneaux de l'ancien Fly's Eye (1981) et la figure 2.10 les tubes photomulti- 
plicateurs (PMT) de Fly's Eye observant la fluorescence U.V. 

Dans cette nouvelle version (HiRes) on a un ensemble de miroirs avec 256 photomul- 
tiplicateurs aux foyers. Chaque miroir observe une portion de 1° du ciel. La lumiere de 
fluorescence se situe dans l'U.V. (300-400 nm). L'utilisation de filtres U.V. permet de 
reduire le bruit de fond. Toutefois la grande sensibilite de chaque cellule est telle que 
l'observation n'est pas possible au clair de lune. La visibility depend de l'energie de la 
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Fig. 2.7 - Detecteur HEGRA. 



gerbe : elle est de l'ordre de 20 km. Elle est limitee par l'absorption de la lumiere par des 
aerosols et des poussieres. Cela donne une surface equivalente a 6000 km 2 . 

L'image de la gerbe qui excite des portions d'atmosphere est enregistree par les cellules 
de photomultiplicateurs. L'energie et la direction de la gerbe peuvent etre reconstruites a 
partir d'un site unique. Toutefois la possibility d'avoir une vue stereo de la gerbe pourrait 
accroitre significativement la resolution en position et en energie. 

2.2.2 Interaction des rayons cosmiques UHE avec la CMBR, la limite GZK 

L'interaction resonante des RC avec les photons du fond a 2.7 K limite le libre parcours 
moyen et rend l'Univers assez opaque, a partir de E GZ k = 10 20 eV. -Eqzk est la limite de 
Greisen Zatsepin Kuzmin pour les protons. 

On obtient cette valeur (exercice) en considerant le seuil de photoproduction de pions, 
dans le processus 

p + 72-7 -> A+ -> N + ir (2.4) 
Pour un angle p-photon egal a 9, on obtient la relation 

E p -pcos9>^^ (2.5) 

ou k est l'energie typique du photon. Pour une temperature de 2.7 K, on obtient 

k = 2.36 10" 4 eV (2.6) 
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Fig. 2.8 - Detecteur HiRes Fly's Eye. 



d'ou la valeur de Eqzk- 

Quelle est la valeur du libre parcours moyen d'un proton d'energie egale a E G zk ? 
Si Ton prend une section efiicace de 10~ 28 cm 2 et une densite de photons de 420 cm -3 , 
on obtient 

Agzk = ~ 8 M P C ( 2 - 7 ) 

Done si a chaque interaction le proton perd environ 20% de son energie, apres quelques 
dizaine de Mpc il se trouvera certainement au-dessous de la limite GZK. On en conclut 
que l'energie des protons a la limite du spectre se degrade tres vite et que leur origine ne 
peut pas etre tres eloignee. 

A titre de comparaison, Virgo, le cluster de galaxies le plus proche de nous, est a 
environ 15 Mpc. 

La probability pour qu'un proton d'energie E GZK traverse 100 Mpc sans collision est 
de 10~ 5 . Finalement, pour qu'il arrive sur terre avec 10 20 eV, le proton devrait avoir en 
moyenne une energie initiale de l'ordre de 100 fois cette valeur. 

2.3 Le modele de Fermi (du "deuxieme ordre") 

Fermi considere en 1949 un modele stochastique dans lequel les particules sont ac- 
celerees par des chocs successifs avec des nuages magnetises (des "miroirs magnetiques") 
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Fig. 2.9 - Tonneaux de l'ancien Fly's Eye (1981). 



comme dans la figure 2.11. 

Considerons un image de masse M et vitesse V et une particule de masse m « M et 
de vitesse v qui tombe sur le nuage sous un angle 9. 

Si (p,E) est le 4-vecteur associe a m, dans le centre de masse du systeme (il coincide 
avec le ref. du nuage car m <<M) le quadrivecteur devient : 

E' = 7(E + Vp x ) 

P' x = l(p x + VE/c 2 ) (2.8) 

avec 7 = (1 — V 2 /c 2 )~ 1//2 et p x = pcosO. 

Apres le choc, l'energie dans le cm. est la meme, seul p x change de signe. Si Ton 
revient au systeme initial : 

E" = 7 (£' - V(-p' x )) = 7 2 [^(1 + V 2 /c 2 ) + 2Vp x ] (2.9) 

Au deuxieme ordre en V/c, on obtient 

V 2 

E" ^E + IE— + 2p x V (2.10) 



E"-E V 2 p x V 2 V 
~^r- « 2— + 2^V = 2— + 2- 

E & E & c 



2— + 2^V = 2— + 2-vcose (2.11) 



25 



UNIVERSITY OF UTAH 




Fig. 2.10 - Les PMT de Fly's Eye. 




FlG. 2.11 - Choc contre un miroir magnetique. 



Cette expression nous permet d'estimer le changement d'energie dans un choc. On va 
effectuer une moyenne sur les angles possibles. On congoit intuitivement que la situation 
d'un choc frontal a plus de chance de se produire, par unite de temps, que celle d'un choc 
depuis l'arriere. Cette probability est classiquement proportionnelle a (v + Vcos^). Si v 
est proche de c, l'expression relativiste est donnee par 

Prob(6) oc 7(1 + cosOV/c) 

Le calcul de la valeur moyenne du terme en cos dans (eq 

/2Vcos6\ _2V J Prob(9)cos6 sin 9d9 
\ c / ~ ~c J Prob(6) sin 6d9 



(2.12) 

. 2.11), nous donne, avec v=c, 



2 ( V x 1 
3 



(2.13) 
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et 

=^ -<¥>"§©' 

Cette expression montre que l'energie varie seulement au deuxieme ordre en V/c. 

Supposons maintenant que la distance typique entre les miroirs soit L, le long de 
lignes de champ magnetique. Le temps entre deux collisions, pour une particule qui suit 
une trajectoire a angle $ par rapport a la ligne, vaut L/(c cos$). La valeur moyenne est 
2L/c, done on estime le taux d'accroissement d'energie moyen a 

dE AE E 8 (VV E 4 (V 2 \ ^ ^ 

M =-rZt = 3UJ Wc = 3\-cL) E = eE (2 ' 15) 

Nous voulons maintenant estimer le spectre en energie E, c. a d. dN/dE. L'equation 
differentielle qui exprime la variation du nombre de particules d'energie E contient les 
contributions suivantes. 



dN(E) 

— = diffusion+gain/perte+sources+fuites (2-16) 

(JjL 

On ne considere pas la diffusion et on suppose qu'on n'a pas de source d'energie E. Le 
terme gain ou perte contient (eq. 2.15), il s'ecrit 



_d_ 

dE 



N(E)§(E) 



(2.17) 



Supposons que le temps de permanence moyen dans la region d'acceleration soit T, inde- 
pendant de E. Les fuites a l'energie E ont typiquement le taux N(E)/T. A l'equilibre ces 
deux contributions doivent s'annuler, pour que la derivee dans (eq. 2.16) soit nulle. On 
obtient : 

ON N(E) ( 1 \ , 

+ ( 2 - 18 ) 



dE E \eT 

ce qui nous conduit, par integration a 

N(E) oc E~ a (2.19) 

avec 

a = 1 + (2.20) 

On a done obtenu une loi de puissance pour le spectre. 
Toutefois un certain nombre de doutes sont licites : 

1) Si Ton considere l'idee des nuages magnetiques on sait que dans notre galaxie les 
vitesses typiques sont de l'ordre de 10~ 4 c. Le parcours moyen dans la galaxie devrait etre 
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d'un pc, ce qui donne une collision/annee. Le gain en energie en un temps raisonnable 
semble hors de portee pour ce modele. 

2) Ce modele ne tient pas compte des pertes d'energie de differentes sortes. 

3) On observe une valeur a ~ 2.5 quelle que soit la direction consideree. Cela est 
difficilement explicable par l'expression (eq. 2.20) qui contient des parametres qui ne sont 
pas des constantes universelles. 

L'essence du processus de Fermi se resume en considerant une probability de fuite par 
collision egale a P. Done, apres k collisions on aurait un nombre de particules 

N = N P k (2.21) 
On considere aussi un changement moyen d'energie par collision 

AE = -i c E (2.22) 
et si Ton part d'une energie E , apres k collisions on aura une energie 

E = E p k avec (3 = 1 + 7 C (2.23) 

On obtient : 

N 
iVo 

Le spectre d'energie est 

N(E)dE = const x E- 1+lnP/ln(3 dE (2.25) 

On voit done que dans le modele de Fermi de deuxieme ordre la valeur de f3 est reliee a 
la variation d'energie par collision et de celle de P au temps T. 



E 

Eq 



lnP/ln/3 



(2.24) 



2.4 Le modele de Fermi (du "premier ordre"), les ondes de choc 

Dans ce modele, l'acceleration des particules cosmiques a lieu dans les regions ou un 
gaz ionise est le siege d'une onde de choc a vitesse supersonique. Cela peut arriver par 
exemple lors de l'explosion d'une Super Nova (SN). 

Les particules qui se trouvent a proximite du front d'onde peuvent traverser le front, 
etre renvoyees en arriere plusieurs fois de suite et gagner globalement de l'energie cinetique. 
Le mecanisme de reflexion sera discute plus loin. On va maintenant se concentrer sur le 
processus d'acceleration et on va tout d'abord resumer quelques resultats de la theorie 
des chocs supersoniques. 

La figure 2.12 montre une onde de choc supersonique qui se deplace dans un gaz defini 
par sa pression p, sa temperature T et sa densite p. Le gaz est au repos devant le front 
d'onde. Si Ton se place dans le referentiel du front, le gaz traverse le front avec la vitesse 
v\ et on le retrouve avec la vitesse V2 a l'arriere. Le traitement du probleme (voir p. ex. 
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gaz au 
repos 



Fig. 2.12 - En haut : une onde de choc supersonique traverse le gaz a vitesse U. Les 
parametres du gaz sont la pression, la temperature et la densite. En bas : on se place dans 
le referentiel du front d'onde. Le gaz traverse le front avec les vitesses indiquees. 



L.D. Landau et E.M. Lifshitz. Traite de physique theorique ; Vol. 6 Hydrodynamique) est 
base sur les ingredients suivants : 

la conservation de la masse impose : 

P1V1 = P2V2 (2.26) 
la conservation de la quantite de mouvement : 

Pi + P\v\ =P2 + P2V2 (2.27) 
Bernoulli implique la conservation de l'energie : 

X -v\ + w x = ^v 2 2 + W2 (2.28) 



sur une ligne de courant, avec 

w = e + pV = t 

P 



w = £ + p V = e- (2.29) 



e est l'energie libre et V le volume specifique. 

Le resultat pour une onde supersonique Mach M=U/ (vitesse du son du gaz au repos) 
dans le cas M » 1 est : 

P2 



Pi 7 + 1 
29 



(2.30) 



P2 _ 7 + 1 

Pi 7-1 



(2.31) 



T2_ _ 2 7 ( 7 -l)M 2 
T\~ (7 + 1) 2 



(2.32) 



avec 7 le rapport des chaleurs specifiques, 7 = 5/3 pour un gaz monoatomique. Dans ce 
cas le rapport des densites vaut 4 et, dans le repere du front d'onde, la vitesse a l'arriere 
est 1/4 de la vitesse a l'avant, cette derniere valant U : 



Si Ton se place dans le repere en aval du front d'onde, ou le gaz est au repos, on verra le 
gaz en amont voyager vers nous a la vitesse V = (3/4) U. 

Une chose a remarquer est que l'epaisseur du front d'onde est en general faible, de 
l'ordre du libre parcours moyen d'une particule du gaz. 

Dans le modele de Bell (1978) et d'autres auteurs, le gaz ionise est le siege d'une onde 
de choc provoquee par une explosion. On suppose que des particules de haute energie 
se trouvant dans le gaz au repos, en aval du front, subissent des diffusions et que leur 
distribution devient ainsi isotrope par rapport au milieu. Dans leur mouvement, certaines 
de ces particules peuvent traverser le front. On se place maintenant dans le repere ou le 
gaz en amont est au repos. On suppose qu'ici aussi la particule subit des diffusions qui 
rendent isotrope les vitesses. La situation est done symetrique, car ici on voit le gaz en 
aval se propager vers nous a la vitesse V = (3/4) U et en traversant le front, la particule 
subit un processus identique a celui decrit plus haut. II n'y a pas de difference entre les 
situations d'un parcours dans le sens amont-aval ou en sens inverse. On peut montrer, par 
des arguments similaires au cas du modele de Fermi de deuxieme ordre, que la variation 
d'energie est donnee par : 



pour une particule relativiste qui fait un parcours complet. Cette variation est done du 
premier ordre en U/c. La raison est qu'ici on n'a jamais de choc qui peut ralentir la 
particule, contrairement au cas vu precedemment. Dans ce modele : 



II nous reste a determiner la probabilite de fuite P. On utilise un argument de Bell. Pour 
les particules de haute energie traversant la surface de choc, le flux par unite de temps 
et de surface est Nc/4, ou N est la densite des particules. En effet, pour une particule 
la probabilite de traverser le front est c x cos9 dans l'intervalle [0, 7r/2] et ailleurs. Le 
nombre de particules dans [9, 9 + d9]d(f> est Nsin9d9d(p et la fraction relative a Tangle 
solide total sin9d9/2 (apres integration sur <f>). Le flux est done : 




(2.33) 




(2.34) 




(2.35) 




(2.36) 



o 
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D'autre part la population des particules de haute energie en amont, apres avoir subi le 
processus de diffusion, sera isotrope par rapport au gaz qui a la vitesse v-i = U/4. Le 
flux perdu sera done NU/4 et la fraction des particules perdues par suite du processus 
d'acceleration elementaire est (NU/4)/(Nc/4)=U/c. On en deduit la probability de fuite 

(? « 1) : 

P = l-U/c (2.37) 
InP = In M - ss — ^ (2.38) 



\n/3 = hi (l + ^j ~ 7 ( 2 - 39 ) 

En utilisant la relation (eq. 2.25) : 

N(E)dE = const x E~ 2 dE (2.40) 

On a done obtenu un resultat "universel", avec un exposant proche de la valeur 2.5 observee 
dans le RC. 



2.5 Le processus d'acceleration, analyse dimensionnelle 

La region des ondes de choc produites dans les supernovae (SN) est une candidate de 
machine acceleratrice. L'image due a Hubble montre la presence d'un front d'ondes de 
choc dans la nebuleuse du Crabe. Cette nebuleuse est nee de l'explosion d'une SN en l'an 
1054. Un calcul dimensionnel montre toutefois qu'on ne peut pas aller au-dela de 1000 
TeV en utilisant les SN. 

Considerons un "accelerateur" de taille R, siege d'un champ magnetique de valeur B. 
Dans les unites appropriees, l'ordre de grandeur du rayon de courbure de la trajectoire 
d'une particule qui a ete acceleree jusqu'a l'energie E vaut : 

R c = E/B (2.41) 

On peut admettre que pour que l'accelerateur fonctionne convenablement, il doit etre 
capable de contenir la particule jusqu'a l'energie maximale ; done il faut que R c <C R '. 

E m ax BR 

10 5 TeV ~ °3 10- 6 G 50pc 

Les unites ont ete choisies de fagon a s'accorder aux valeurs attendues dans des chocs 
associes aux SN. 

On peut supposer que le genou a 10 15 eV est associe au "cutoff" abrupt de la compo- 
sante des rayons cosmiques generee par les SN. Pour aller au dela, on doit immaginer soit 
des regions acceleratrices beaucoup plus grandes soit des champs plus intenses. 

Les candidats les plus probables a l'acceleration sont les noyaux de galaxies actives 
(AGN). 
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(2.42) 
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Fig. 2.13 - Nebuleuse du Crabe. 



Dans des cas extremes, on considere que les etoiles a neutrons et les trous noirs en 
collision sont responsables d'episodes explosifs avec une energie suffisante. Ces evenements 
seraient associes a la production de "gamma-ray burst" (voir plus loin). 

A-t-on des preuves de l'exactitude du modele de Fermi ? Un exemple suggestif, bien 
que les ordres de grandeurs soient completement differents, est l'acceleration des particules 
dans le vent solaire. A l'occasion d'une onde de choc observee par un satellite, on a mesure 
des particules energetiques un peu avant et un peu apres le choc. Les resultats obtenus 
indiquent effectivement la presence d'un processus d'acceleration par onde de choc, comme 
decrit plus haut. 

2.6 Les "flammes" du Soleil 

Ce sont les evenements les plus energetiques observes dans notre Soleil. La luminosite 
du Soleil, dans la ligne Ha en particulier, augmente enormement sur une surface de 
10 9 km 2 . Lors des eclipses, on observe des filaments qui traversent la chromosphere et 
qui s'elevent tres haut sur la surface. Le phenomene dure environ 1 heure avec souvent 
un flash initial, la phase impulsive, d'une duree de quelques secondes, observable dans 
les micro-ondes, dans les U.V. et X durs. Cette impulsion est attribuee a l'acceleration 
d'electrons, sur une echelle de temps tres courte. 

L'energie liberee est de 10 22 a 3 10 25 J. Le diametre est L=3 10 7 m, la hauteur a ete 
estimee a H=2 10 7 m, ce qui nous donne un volume de 2 1 22 m 3 . L'energie thermique 
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contenue dans ce volume n'est certainement pas suffisante pour expliquer le phenomene. 
Par contre. si Ton considere le contenu magnetique, l'energie totale dans le volume V 
est E = V(B 2 /2fx ) = 10 28 B 2 J; pour B=0.03 a 0.1 Tesla, on obtient le bon ordre de 
grandeur. Ces valeurs de B sont celles observees lors de remission micro-onde dans la 
phase impulsive. 

L'etude de ce phenomene nous montre un cas interessant de processus de rechauffement 
du plasma et un modele possible d'acceleration de particules. L'etude de la relation entre 
le mouvement du plasma et la ligne de champ B montre que ce champ suit la matiere. 
Le champ B est condamne a suivre le mouvement du plasma, qui peut etre capricieux. 
Par exemple, le mouvement convectif peut amener le champ B a prendre localement des 
valeurs extremement elevees. 

On rencontre des regions ou le champ B a sa direction inversee. C'est dans ces regions 
que l'energie magnetique a la possibility de se dissiper, si le mouvement de convection 
rapproche les lignes de champ et si le plasma dans la region neutre a une resistivite r/ 
suffisante. Cette derniere peut se calculer a partir de la theorie des plasma 

7] ss 10 3 /T 3/2 ohm m (2.43) 

Dans la figure 2.14 on montre une region ou le champ B est dirige selon l'axe x et sa 
polarite est inversee entre les regions y>0 et y<0. Dans ce cas une densite de courant 
elevee parallele a z est presente. On suppose le systeme infiniment etendu selon l'axe z. 




Fig. 2.14 - Les lignes du champ B sont dirigees en sens oppose. Elles s'approchent, et 
vont finalement se reconnecter. 



Le plan xz constitue done une bande magnetiquement neutre. Supposons que les lignes de 
champ se rapprochent en suivant le mouvement de convection a vitesse v, v//y, comme 
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indique dans la figure 2.14. Pour preserver la continuity, le plasma et le champ B qui 
se trouvent dans cette bande sont obliges de sortir le long de Taxe x. Puisque le plasma 
possede une resistivite ohmique 77, ce mouvement entraine une dissipation d'energie. 
Considerons une boucle qui entoure la region et appliquons le theoreme de Stokes : 

JdS = — lB-dl (2.44) 
Mo J 

Si t est la largeur de la boucle et L sa longueur, on obtient les ordres de grandeur : 

iLJ « — 2BL J « — (2.45) 
pL Hot 

(ce qui correspond a J = rotB//x, dans la region interessante) . 
L'ordre de grandeur de la puissance dissipee par unite de volume est : 

" j2k "(^) 2 (2 - 46 > 

Avec la diminution de la largeur £, la densite de courant augmente quadratiquement et 
le processus devient de plus en plus efficace. Comme valeur minimale de £, nous pouvons 
prendre le rayon de courbure des particules dans B. 

Dans ce calcul, nous devons aussi tenir compte de la valeur de la pression. Celle- 
ci est definie par l'egalite entre la pression du gaz et celle du champ B. A l'equilibre la 
compression est compensee par l'echappement de matiere selon x. Si l'on neglige les forces 
magnetiques, le mouvement du plasma selon x est decrit par : 

dp dv x ( d _ „\ , 

-i=^= s "U + "' v )^ (2 - 47) 

A l'etat stationnaire dv x /dt = 0, 

dp dv x 



(2.48) 



dx dx 
Par integration entre x=0 et l'infini, on obtient : 

Ap = X -pvl (2.49) 
Cette difference de pression doit etre compensee par celle associee au champ B. 

Ap = |f (2.50) 

Cela nous permet de calculer la vitesse de fuite du plasma selon x : 

v x - = v A (2.51) 
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v A est la vitesse d'Alfven analogue a la vitesse du son dans un milieu magnetohydrody- 
namique. Par continuity, la vitesse d'injection vaut : 




v = v A - (2.52) 

Si Ton considere que seul Tapport d'energie sous forme magnetique est important, le 
flux d'energie vers la region de dissipation vaut ~ vB 2 . On egalise Tapport d'energie 
magnetique au taux de dissipation ohmique : 

f f B 2 

v J 2 dV = J J T^vdS (2.53) 

Sur une longueur unitaire en z on obtient : 

B 2 B 2 

r]J 2 £L = 2 vL 7]J 2 i = —v (2.54) 

2/x /i 

Le facteur 2 correspond a la surface d'en haut plus celle d'en bas. En combinant avec (eq. 
2.46) 

v = 477/^0 (2.55) 
En utilisant (eq. 2.52), on peut resoudre par rapport a £ et v : 

t 2 = (2.56) 

/J, V A 

v 2 = -%v A (2.57) 
fi L 

L'echelle de temps de la reconnexion des lignes est : 

t R w i/v = L/v A (2.58) 

Dans le cas du soleil, juste avant une flamme : 

T = 2 10 6 K (2.59) 

L = 10 7 m 

B = 0.03 T 

n = 10 16 m- 3 

v A = 6 10 6 m/s 

t R pa 1 s et la vitesse de reconnexion vr ~ 0(1) m/s 

L'energie dans la bande neutre est de 3 10 16 Joules, la luminosite est done 3 10 16 W. Si 
Ton prend l'energie totale d'une flamme petite/moyenne, 10 23 J, on deduit un temps de 
dissipation de l'ordre de 10 7 s. Cela est en complet desaccord avec l'observation, car ces 
phenomenes ont lieu sur une echelle de temps de l'ordre de l'heure. 
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Pour des modeles alternatifs voir la ref. 1) mentionnee au debut du cours. 

Lors de l'apparition des flammes solaires on a observe des protons de plus de 200 MeV, 
et des electrons de 10 MeV. 

Dans la couche de reconnexion, les valeurs de rotB et dB/dt sont certainement non 
nulles, ce qui induit un champ E accelerateur. La question est d'estimer l'efficacite d'un 
tel modele d'acceleration. Si Ton prend les ordres de grandeurs : 

V x E = -dB/dt =>- E/l « B/t R E « Bv R (2.60) 

on obtient un champ E de l'ordre de 0.01-0.1 V/m. Sur une distance de 10 7 m, on peut 
envisager des accelerations de 1 MeV au maximum. 

II apparait que meme dans le cas des flammes solaires des processus "a la Fermi" 
doivent etre envisages, si Ton veut rendre compte des energies mesurees. 



2.7 Acceleration dans la magneto-sphere des pulsars 

Ce modele d'acceleration est peut etre la seule alternative au modele de Fermi. Les 
pulsars possedent un champ magnetique tres fort, jusqu'a 3 10 8 T: ils tournent tres vite, 
leur periode varie entre la ms et la dizaine de secondes. Une certaine quantite d'energie 
est dispersee dans l'espace, ce qui provoque un lent mais inexorable ralentissement de la 
rotation. 

Dans ce genre de modele, on considere en general un aimant en rotation, avec l'axe 
magnetique non aligne sur la direction du champ. Toutefois, si Ton considere le cas simple 
d'un alignement exact, on peut apprecier l'importance des effets e.m. par rapport a la 
gravitation. 



F Pm ev x B ettBr 3 

(2.61) 



em 

r^j 



F grav GMm/r 2 GMm 

Si Ton prend le cas du Crabe, par exemple, la periode vaut P=27r/f2=33 ms. Pour un 
electron en orbite autour du pulsar, le calcul donne un ordre de grandeur 10 12 , ce qui 
montre que tout ce qui se passe a la surface du pulsar est domine par les forces e.m.. 
La d.d.p. entre le pole et l'equateur d'un pulsar avec P = 1 s, B = 10 8 T est de 10 16 
V. Les electrons sont done arraches de la surface et rapidement acceleres. Par suite des 
interactions, ces particules vont generer une veritable gerbe e.m.. C'est probablement le 
phenomene qui produit la luminosite des pulsars. Par exemple, toujours dans le Crabe, on 
observe une luminosite equivalente a 10 11 K, ce qui correspond a des electrons d'energie 
d'environ 10 MeV. 

De fagon generale, des arguments dimensionnels montrent qu'il est possible que les 
etoiles a neutrons soient responsables de l'acceleration des rayons cosmiques, jusqu'aux 
energies extremes. L'ordre de grandeur de la limite superieure de l'energie peut s'obtenir 
en considerant une etoile de dimension L avec un champ B. Les variations de flux peuvent 
avoir lieu sur une echelle de temps minimale 5t = L/c. On en tire la relation : 

Vxl= -dB/dt E/L^B/5t E ss Be (2.62) 
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L'energie acquise par la particule donnee par eEL = ecBL. Si Ton prend des valeurs 
appropriees pour la surface d'une etoile a neutrons L=100 km, B=10 6 T, on obtient une 
energie de 10 19 eV. Done une etoile a neutrons a a sa disposition l'energie necessaire pour 
accelerer les particules aux energies extremes. II faut toutefois elaborer ce modele pour 
prouver que cela est effectivement possible. 

2.8 Au dela du Modele Standard 

Au-dela du processus d'acceleration, on peut imaginer que la portion extreme des RC 
soit d'origine cosmologique et associee a la disintegration d'objets exotiques. 



i — i — t — i — | — i 'f t i — | — i — i — i — i — ] — r -it "I i 




LOG (Energy [eV]) 

Fig. 2.15 - Prediction theorique pour le spectre des cosmiques UHE a partir du modele 
des colliers (courbes en trait continu). 



Dans le modele du Big Bang, l'Univers traverse plusieurs transitions de phase, associees 
a des changements de symetrie. Par exemple, au cours du processus de refroidissement, 
au moment ou la temperature etait de l'ordre du TeV, le mecanisme de Higgs intervient 
pour donner une masse aux bosons intermediaires W et Z, tout en laissant le photon sans 
masse. G'est le champs de Higgs qui reduit la symetrie de l'Univers selon : 

SU(2) x U(l) -> U(l) em (2.63) 

La transition de phase se fait parce que le champ de Higgs prend une orientation bien 
determinee. Toutefois, des regions de l'espace causalement deconnectees pourraient etre 
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le siege de champs orientes differemment. II s'agit d'une situation analogue p. ex. aux 
domaines magnetiques qui se forment au moment ou une substance ferromagnetique est 
refroidie en dessous de sa temperature critique. 

Plusieurs sortes d'objets sont envisageables, avec beaucoup d'energie stockee. A une 
dimension, on parle de monopoles. Les cordes sont a 2D, les domaines a 3D et les textures a 
4D. Des topologies plus compliquees ont ete considerees, comme les "colliers", qui donnent 
les predictions theoriques de la figure 2.15. 
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3 La radiation e.m 



3.1 Introduction 

La premiere source d'information sur le Cosmos est la lumiere observee directement ou 
par des telescopes. Les instruments optiques les plus performants en 1999 sont le telescope 
de Hubble (en orbite) et le Very Large Telescope (VLT, observatoire du "Mont Palomar" 
au Chili). 

L'observation de la radiation de fond a 2.7 K se fait a l'aide de radiotelescopes. La 
longueur d'onde typique de ce spectre est de quelques mm. 

Le signal des pulsars a ete decouvert dans le domaine radio. II en est de meme des 
quasars. Une certaine classe de Galaxies est tres active dans cette region du spectre. 

On a mis en orbite des satellites avec des detecteurs X et gamma. 

Globalement, le domaine e.m. d'observation couvre plus de 21 ordres de grandeur 
d'energie (de 10~ 5 eV a 10 15 eV). La figure 3.1 montre ce domaine, ainsi que les moyens 
experimentaux a disposition de l'observateur. Dans le domaine 100 MeV<E<10 GeV, on 
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Fig. 3.1 - Domaine electromagnetique sonde par les astrophysiciens. On a indique les 
noms des dispositifs "type" utilises dans la region gamma. 



etudie la radiation gamma par des satellites. A plus grande energie, seuls les detecteurs au 
sol ont une acceptance et resolution suffisantes. L'interaction avec l'atmosphere produit 
une gerbe e.m. qui peut etre detectee par l'observation de la fluorescence ou de la lumiere 
Cherenkov. La phenomenologie est tres riche : on peut observer des spectres continus 
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ou des raies en emission ou absorption. Les emissions sont continues, pulsees (pulsars) ou 
explosives (Gamma Ray Bursts). 



Annihilation 
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Fig. 3.2 - Mecanismes de production de photons. 
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3.2 Mecanismes de production 



La figure 3.2 montre les mecanismes principaux de production de photons. 

Dans l'etude du rayonnement, il faut tenir compte des effets qui modifient le spectre, 
en particulier l'effet Doppler du au mouvement de la source, le "red shift" cosmologique 
et le "red shift" gravitationnel. Ce dernier vaut : 

GM E-y , . 

Par exemple, a proximite d'une etoile a neutrons, on peut avoir facilement un "red shift" 
de 10-30%. 



3.2.1 Les raies spectrales 

La figure 3.3 montre un spectre obtenu en regardant en direction du centre Galactique. 
Les raies spectrales dans la plage d'energie des gamma peuvent etre produites par les 
phenomenes suivants : 




Fig. 3.3 - Spectre gamma obtenu en regardant en direction du centre Galactique. 



a) L'annihilation e + e . 

Si les deux particules sont produites presque au repos, elles forment un systeme dit 
"positronium" dans un etat 1 S'o pour 25% et dans un etat 3 Si pour les 75% restants. Le 
premier etat s'annihile en deux photons de 0.511 MeV, le second en trois photons dont 
l'energie est distribute dans un continuum de valeur maximale 0.511 MeV. L'annihilation 
en vol, qui sera traitee plus loin, contribue aussi au continuum. Les positons proviennent, 
par exemple, de la chaine pion — > muon — > positon. Dans la nucleosynthese on a aussi 
une production de positons (voir Table 3.1) : 
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57 Co — > 57 Fe donne des e + d'environ 100 keV 

26 Al -> 26 Mg donne des e + d'environ 1 MeV 

Tab. 3.1 Examples de production de positons dans nucleosynthese 



La figure 3.4 montre le spectre observe au voisinage de 511 keV. La largeur de la raie 
est superieure a la resolution experimental, ce qui pourrait indiquer une production en 
vol. 



ENERGY (keV) 




Fig. 3.4 - Raie a 511 keV du centre Galactique (satellite HEAO-3). Une fois corrigee 
par la resolution de l'appareillage, la largeur intrinseque de la raie est d'environ 2 keV, a 
mi-hauteur. 



b) La capture radiative de neutrons. 

Elle produit un noyau compose qui se desintegre par emission gamma. 

Par exemple n + p — > 7+ 2 H produit une raie a 2.2 MeV. La section efficace de ce 
processus vaut 7.3 10~ 20 t> 1 cm 2 (v en cm/s). Le processus ne se realise de fagon efficace 
qu'en presence de neutrons a faible vitesse (neutrons thermiques). 

Un autre exemple est la capture neutronique par le 56 Fe qui produit du 57 Fe ; celui-ci 
se desexcite par emission de gamma de 7.632 et 7.646 MeV. 

Les neutrons sont produits par des reactions du type 4 He(p,pn) 3 He, 3 H(d,n) 4 He. Le 
neutron ayant un temps de vie de 15 minutes, le processus radiatif ne peut se realiser que 
dans des regions a densite de l'ordre de 10 16 protons/cm 3 , ou le fibre parcours moyen est 
plus petit que l'echelle spatiale associee au temps de vie (flammes, disques d'accretion,...). 
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c) Les raies de desexcitation nucleaire . 



16 2.741, 6.129, 6.917, 7.117 MeV 

12 C 4.438 MeV 

M N 2.313 et 5.105 MeV 

16 6.129 MeV etc... 

Tab. 3.2 - Quelques raies de desexcitation nucleaire. 



d) Les niveaux de Landau . 

En presence d'un champ magnetique tres intense, des raies d'emission peuvent etre 
produites par des transitions entre des niveaux de Landau (ou cyclotroniques). Pour le 
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Fig. 3.5 - Raies d'absorption de GRB 880205 observees par GINGA (satellite astrono- 
mique japonnais), attribuees a des niveaux de Landau. 



premier niveau l'energie vaut : 

E~ = h—B (3.2) 

m 

ou ^ = charge electrique/masse de la particule. 

Pour l'electron -S- = 1.76 10 11 fs -1 T -1 l et on obtient : 

£ 7 [keV] = 11.6 x (^p) = H.6 x B 8 [T] (3.3) 

Des raies d'absorption ont ete observees dans Her X-l, dans les GRB, ... Un exemple est 

donne dans la figure 3.5. Pour des energies dans l'intervalle 20-50 keV, le champ B est 
done de l'ordre de 10 8 T. 
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3.2.2 Les continuum 

La figure 3.6 montre un spectre obtenu en regardant en direction du centre Galactique. 




ENERGY (MgV) 

Fig. 3.6 - Spectre en direction du Centre Galactique, observe par COS B et SAS II. 
Les courbes en trait continu donnent les predictions theoriques pour les contributions de 
la radiation provenant du Bremsstrahlung, de Compton et de n°. La courbe en traitille 
donne la contribution totale. 



Les sources de gamma qui contribuent a la partie continue du spectre sont : 

a) I' annihilation en vol d'un positon et d'un electron (un des deux au repos, normalement 
l'electron). La section efficace a grande energie est toutefois assez faible : 

o = ^ [ln(2 7 ) - 1] (3.4) 

avec r e = e 2 /m e c 2 le rayon classique de l'electron et 7 le facteur de Lorentz. A 100 MeV on 
obtient a = 6 10~ 27 cm 2 , et un libre parcours moyen de 1.6 1 26 /n cm, ou n est la densite 
electronique. Pour le milieu interstellaire n est de l'ordre de 3 10 ~ 2 cm -3 , ce qui donne 
environ 2 Gpc. Done cette contribution pour etre efficace necessite des regions denses, 
mais pas trop, car le risque d'autoabsorption augmente aussi avec la densite. 
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Fig. 3.7 - Taux de perte d'energie en fonction de l'energie des electrons. Plusieurs pro- 
cessus contribuent a la production de X ou de gamma. S : emission synchrotron dans un 
champ B=4 /zG ; C„ : effet Compton inverse avec des photons dans le visible dont la den- 
site d'energie vaut 0.45 eV/cm 3 ; C3, idem avec les photons CMBR; B : bremmstrahlung 
pour une densite de 1 atome/cm 3 ; I perte par ionisation ; E : "fuite" avec un temps de vie 
moyen dans la Galaxie de 3.3 10 6 ans. 



b) la radiation de Bremsstrahlung dans le champ coulombien des noyaux et d'autres elec- 
trons. Dans le cas des noyaux, on trouve (B. Rossi, High-energy particles 1952) : 



dE 

a{E e , E 1 )dE 1 = 4a em Z 2 r 2 e —tE(E e , v 



(3.5) 



ou v = E y /E e . 

F est une fonction qui depend de l'environnement. Dans le cas de noyaux completement 
ionises (protons, par exemple), on a : 



F(E e ,u) 



l + (l-z/) 2 --(l-z/) 



In 



f2E e l-v 1 
\mc 2 v 2 



(3.6) 



La forme de cette contribution au spectre est done approximativement proportionnelle 
a E" 1 . On peut exprimer la perte d'energie des electrons par ce processus en considerant 
la longueur de radiation du milieu, X . Sur une longueur egale a X , l'energie de l'electron 
est reduite approximativement d'un facteur 1/e. 
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Fig. 3.8 - A gauche : la particule emet de la radiation synchrotron de preference vers 
l'avant. A droite : l'observateur voit l'intensite de la radiation augmenter quand la parti- 
cule vient a sa rencontre. 




Fig. 3.9 - A gauche : structure temporelle de remission. A droite : spectre des frequences. 



Pour l'hydrogene X =61 g/cm 2 . Dans le milieu interstellaire on a un proton/cm 3 , ce 
qui nous donne Xo=10 Mpc. Done pour v=c, on obtient un libre parcours moyen de 4 10 7 
annees lumieres. La figure 3.7 donne le taux de perte d'energie en fonction de l'energie 
des electrons pour le milieu interstellaire. 

On trouve que l'emissivite due au bremmsstrahlung dans le gaz interstellaire que Ton 
trouve dans le plan galactique, pour un seuil de 100 MeV, est d'environ 2 10~ 25 photons 
cm s . 

Q. : Montrer que, si Ton considere un parcours de l'ordre de 10 kpc, le taux mesure 
par un detecteur qui a un seuil a 100 MeV, une surface de 500 cm 2 , un angle solide de 
0.25 sr et une efficacite de 10% est d'environ 500 gamma/jour. 

c) L'effet Compton inverse. 

Considerons un photon de faible energie, e=0(leV), qui est frappe par un electron 
relativiste. Du point de vue de ce dernier, le photon a une energie et lors de la diffusion 
Compton inverse il aura une energie plus petite que cette valeur. Done, si Ton retourne 
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au systeme de l'observateur initial, on aura un photon d'energie 



E 1 pa e^y 2 quand 75 <C m e c 2 (3.7) 



E 1 pa E e pour 75 > m e c 2 (3.8) 

Par exemple, si Ton considere les photons de 6 10~ 4 eV du CMBR, la transition entre 
les deux regimes se fait a environ 5 10 14 eV. Si Ton considere un spectre d'energie des 
electrons de la forme E~ a , alors l'emissivite due au bremsstrahlung suit une loi en : 

N(E 7 ) oc E^ a+1 ^ (3.9) 

d) L'emission synchrotron dans un champ magnetique. 

L'electron qui spirale dans un champ magnetique emet un rayonnement, de preference 
en avant, comme dans la figure 3.8. La theorie e.m. nous dit que la largeur du lobe est de 
l'ordre de : 

0pamc 2 /E = I/7 (3.10) 

On peut imaginer que les photons sont emis a l'interieur d' un cone d'ouverture 6 autour 
de la trajectoire de la particule. 

L'observateur place a distance considerable voit un champ e.m. arriver par bouffees, 
a la frequence de Larmor (voir la figure 3.9). 



eB 28 



Vn 



B [GHz] avec B exprime en Tesla (3. 11) 



9 2nmy 7 
La duree d'une bouffee vaut environ : 

Q 

At pa — T x (1 - v/c) (3.12) 

ou T est la periode du mouvement cyclotronique de la particule et v la vitesse de cette 
particule. 

Le terme entre parenthese tient compte du fait que l'electron et les photons observes 
se propagent dans la meme direction . Le spectre qui en resulte est riche en harmoniques, 
avec une valeur maximale de l'ordre de 1/T. La valeur de l'energie qui correspond au 
maximum d'intensite est donnee par : 

E max (eV) = 5 10- 13 5 ±7 2 (3.13) 

B± = Bsma en Tesla, a est Tangle entre la vitesse de la particule et le champ. 

L'energie des photons attendue pour ce processus est done assez faible, sauf si les 
champs et les energies sont tres eleves. L'observation de remission de synchrotron indique 
done la presence d'electrons relativistes. Le spectre que l'on obtient, pour une distribution 
d'energie des electrons en E~ a , est exactement comme pour l'effet Compton inverse : 

N(E y ) oc Eltt a+1 W (3.14) 
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La lumiere emise par un electron individuel est polarisee circulairement, mais celle pro- 
duite par une distribution isotrope de particules donne une resultante polarisee lineaire- 
ment, dans le plan orthogonal a la direction du champ magnetique. 

e) La radiation due a la courbure de la trajectoire dans la magnetosphere des pulsars. 

A proximite d'un pulsar, le champ magnetique intense oblige les particules a voya- 
ger selon les lignes de champ. Puisque ces dernieres sont courbees, on obtient un lobe 
d'emission en avant. L'energie typique de la radiation vaut : 

T r \ 3hc-f 3 3 10~ 5 , 
EieV) « - — - = — — -7 3.15 
v 1 2 p p(cm) v 1 

p est le rayon de courbure de la trajectoire dans le champ magnetique. Pour un electron 
de 10 13 eV, un rayon de courbure de 10 8 cm, on obtient E=2.5 GeV. 



Cas de protons. 

Les memes mecanismes que ceux discutes ci-dessus existent pour les protons. Les 
processus de b) a e) sont en general moins efficaces qu'avec les electrons, a cause de la 
masse du nucleon qui est environ 2000 fois plus grande. 

f) Interaction des protons avec la matiere et production de pions neutres. 

La production de pions neutres peut se faire par collision de protons avec la matiere 
interstellaire (protons). Le seuil pour la production d'une particule de masse m (exercice) 
est : 



T seuil = 2mc 2 (l + m/ AM) (3.16) 

ou T est l'energie cinetique du proton et M sa masse. 

Pour la production d'un n° le calcul donne T seui i 280 MeV. 

Le pion neutre se desintegre presque instantanement (0.83 10~ 16 s) en donnant deux 
gamma, 7r° — > 77. Dans le referentiel du 7r° l'energie se repartit egalement entre les deux 
photons. Le boost de Lorentz du pion determine l'energie observee dans le laboratoire. 
Considerons une famille de pions se deplagant a une vitesse v par rapport a l'observateur. 
Suivant l'angle du photon par rapport a la vitesse du pion, l'observateur va mesurer une 
energie affectee par le boost de Lorentz. 

On trouve (exercice) que la distribution est plate, avec les extremes donnes par : 



mc 2 1-3 ^ mc 2 1 + 3 . 



La figure 3.6 montre la contribution de cette composante, observee en direction du 
Centre Galactique. 
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g) Annihilation proton-antiproton p + p — > n n x 71 

Pour une annihilation au repos (protonium) les lois de conservation interdisent n,,- <3. 
La conservation de l'energie implique que n n < 14. On peut montrer (exercice) que pour 
n n =3 l'energie maximale d'un pion est : 

E nma x = (4:M 2 - 3m 2 ) /AM « 923 MeV (3.18) 

Si Ton considere des pions neutres on aura, d'apres (eq. 3.17), des energies de photons 
comprises entre 5 et 919 MeV. La figure 3.10 montre cette distribution. 
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Fig. 3.10 - Spectre d'energie des gamma issus de l'annihilation au repos de protons et 
d'antiprotons. 



h) La radiation du corps noir. 



Cette radiation est presente des que le corps a une profondeur optique elevee pour 
toutes les frequences. Le spectre est donne par la relation de Planck : 

2vr/iz/ 3 " 



I{v)du = 



1 



exp{his / kT} — 1 



du MeV cm s 



2 „-i 



N{E 1 )dE 1 = 



I{v)dv 
hv 

9.89 10 40 



(3.19) 
(3.20) 



exp{1.16 W 10 hu/T} - 1 



2 „-i 



E^dE 7 photons cm s 



La loi de Wien donne l'energie au maximum de la distribution : 

E wieQ = 4.7 10" 10 T ou T est en degre K et E en MeV 



(3.21) 
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L'energie moyenne des photons est : 



E = 2.7kT = 2.3 10 



-10 



T 



(3.22) 



Pour obtenir des photons de l'ordre du MeV, il faut done un systeme a 10 10 K. 
Q. : Estimer la fraction d'energie de l'Univers, f2, que l'on trouve sous forme e.m.. 

3.3 Mecanismes d'absorption 

On a essentiellement deux cas : l'interaction des gamma avec la matiere et l'interac- 
tion des gamma avec un champ e.m.. Ce dernier peut etre un photon reel ou un champ 
magnetique. 

a) Interaction gamma-matiere : il s'agit de l'effet Compton, de l'effet photoelectrique et 
de la creation de paires. 

L'effet Compton et l'effet photoelectrique ont une importance significative en dessous 
de 10 MeV pour de la matiere a faible Z, ou en dessous de 100 MeV pour un Z eleve. En 
dessus de ces energies, la creation de paires e + e~ est dominante. 

bl) Interaction gamma-gamma. 



Pour un photon d'energie £i qui frappe un photon cible d'energie e 2 sous Tangle 
d'incidence 9, l'energie totale disponible pour la creation d'une paire de particules de 
masse m, est donnee par : 




Target Photon 



E= ^2 £l e 2 ^ 



cos 9) 



(3.23) 



done le seuil est a : 



2m 2 c 4 



(3.24) 



£2(1 — cos#) 



Dans le cas d'un choc frontal, le seuil de production d'une paire e + e vaut : 



£1 = 0.26 I0 12 /e 2 



ou les energies sont exprimees en eV 



(3.25) 
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On peut calculer (exercice) les valeurs du seuil de production e + e pour les cibles 
canoniques : 



CIBLE E seuil 

2.7 K 400 TeV 

Starlight 0.1 TeV 

1 keV X 0.3 TeV 

Tab. 3.3 - 



A tres grande energie, (ou = ^/£iS2 3> m), la section efficace vaut 

2 



2 f m \ 



21n|^-l 

m J 



avec r e = e 2 /mc 2 = 2.82 x 10 15 m le rayon classique de l'electron . 
Quand m est grand par rapport a u 



a = nri [ 1 - -f 



(3.26) 



(3.27) 
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Fig. 3.11 Interaction de gamma avec des cibles photoniques. BB est la radiation CMBR. 
La creation de paire de muons se manifeste vers 10 20 eV. 



L'efficacite d'un tel processus depend de la densite de la cible. Dans le cas des photons 
CMBR (400 par cm 3 environ) la longueur d'attenuation d'un photon de 10 15 eV est de 
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l'ordre de 10 23 cm. La figure 3.11 montre la longueur d'attenuation pour differentes cibles, 
en fonction de l'energie du photon. 

Si l'on considere une distance equivalents au rayon de Hubble, on peut dire que l'Uni- 
vers ne devient transparent qu'en dessous de 10 12 eV. 

b2) Interaction gamma-champ magnetique. 



En presence d'un champ magnetique intense, le seuil de production e + e est donne 
par la relation suivante : 



Done un photon de 40 MeV peut se convertir dans le champ d'une etoile a neutrons de 
B=10 7 Tesla. 

3.4 Les detecteurs 

Les caracteristiques en sont : le domaine d'energie, la resolution en energie, l'acceptance 
et la resolution angulaires, l'efficacite et la surface de detection. 

La resolution temporelle est aussi un atout pour pouvoir resoudre la structure tempo- 
relle d'un episode Gamma Ray Burst (GRB) et/ou pour la rejection du bruit. Du point 
de vue du mecanisme de la detection, on a 3 domaines d'energie dans l'interaction gamma 
matiere : jusqu'a 3 MeV l'effet photoelectrique domine ; de 3 a 10 MeV l'effet Compton 
domine ; la creation de paires e + e~ est dominante au-dela de 10 MeV. 

3.4.1 Les detecteurs dans l'espace 

II s'agit de detecteurs places sur des ballons, des fusees, des navettes spatiales... L'ins- 
trument de la NASA, appele "Compton Gamma Ray Observatory" (CGRO, figure 3.12), 
en service de 1991 a 2000, en est un exemple. II comprenait 4 detecteurs specialises cou- 
vrant le spectre electromagnetique de 30 keV a 30 GeV. Dans l'ordre croissant de la plage 
d'energie couverte, les detecteurs etaient : 

1) le "Burst And Transient Source Experiment" (BATSE) avec 2 modules comprenant 
2 scintillateurs Nal(Tl) chacun places aux 4 coins du satellite, examinant la totalite du 
ciel ; 




e 



-^seuii(eV) = 4 10 W /B avec B exprime en Tesla 



(3.28) 
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Fig. 3.12 - Les detecteurs du "Compton Gamma Ray Observatory". 



2) 1' "Oriented Scintillation Spectroscopy Experiment" (OSSE), avec des scintillateurs 
(NaI(Tl)+CsI(Na)) equipes de collimateurs ; 

3) le "Compton Telescope" (COMPTEL) utilise Teffet Compton dans des scintillateurs; 

4) 1' "Energetic Gamma- Ray Experiment" (EGRET), un detecteur avec chambres a 
etincelles, pour des gamma de 20 MeV a 30 GeV. 

• A basse energie (< 3 MeV) on utilise un scintillateur couple a un photomultiplica- 
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teur. Pour attenuer les effets de la radiation parasite, on doit prevoir un ecran. Par exemple 
on peut entourer le tout avec un ecran en plomb qui definit aussi Tangle solide observe par 
le scintillateur. Au prix d'une reduction de l'efficacite, on peut prevoir un ecran "actif" 
(figure 3.13) qui est place en "anticoincidence" avec le detecteur. 

ACTIVE 




CENTRAL 
DETECTOR 



Fig. 3.13 - Telescope gamma. On utilise un cristal Nal couple a un PM. La "collimation" 
se fait par un ensemble de cristaux qui definissent une ouverture de 15°. 




Fig. 3.14 - Diffusion Compton. 



• Pour le domaine 3-10 MeV l'effet Compton peut etre utilise pour etablir la direction 
de la source, sur une base statistique. 

La figure 3.14 montre un photon d'energie Eo qui entre en collision avec un electron 
cible a T arret (masse m). Le photon est diffuse sous Tangle 6 avec une energie Ei. On 
trouve : 

Ei = r (3.29) 
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Fig. 3.15 - do jd9 et da/d(f> pour la diffusion Compton du photon et rejection de l'electron 
de recul. 



L'electron de recul est ejecte avec une energie E telle que : 



= 4 L+£ggg (3.30) 



La relation entre les angles d'emission est : 



cot0 = I 1 + ) tan(0/2) (3.31) 



mc 2 

La section efficace differentielle (formule de Klein Nishina) pour la diffusion d'un photon 
non polarise est : 

^) = _f!_f|iyf^ + |l_ sin2S ) ,3.32) 
dVt 2mc 2 \Eq) \E l E J y ' 

Pour un photon emis dans la fenetre [9, 9 + d9], l'electron est emis dans [0, <f> + d<f>]. On 
peut done transposer l'equation 3.32 dans une section efficace exprimee pour l'electron de 
recul. 

Les deux figures 3.15 montrent da/d9 a gauche et da/dtp a droite pour le photon et 
l'electron respectivement, pour differentes energies du photon incident (a = E /mc 2 ). 

Un telescope Compton peut utiliser le photon diffuse ou l'electron de recul pour donner 
une estimation de la direction du photon entrant. Le systeme de base qui utilise le photon 
sortant est donne a la figure 3.16. A est un scintillateur plastique qui est place en anti- 
coincidence dans la procedure de declenchement. Cela permet la rejection des cosmiques 
charges. Le photon diffuse dans S est observe par C. A e est aussi en anticoincidence et 
permet d'eliminer les evenements avec un electron diffuse vers C. Le signal de "trigger" 
est done : tviggev=AA e SC . D'apres les relations vues plus haut, on peut estimer Tangle 
9, si Ton connait les energies en jeu : 

Eq — El n E n 

cose* = 1 — — — mc =1 — W\ mc (3.33) 

E\Eq EiyE -\- EiJ 
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Fig. 3.16 - Configuration de base pour un Fig. 3.17 - Telescope du Max Planck Insti- 

telescope Compton qui se sert du photon dif- tut de Munich, 
fuse pour donner une estimation de la direc- 
tion probable du photon incident. 

II faut done connaitre l'energie du photon diffuse. Idealement C est un calorimetre. 
Un exemple est represente a la figure 3.17. II s'agit du telescope du Max Planck Institut 
de Munich. COMPTEL a essentiellement la meme configuration. Les scintillateurs d'an- 
ticoi'neidence ne sont pas montres. Chaque evenement i definit un cone de demi-angle 0j. 
Pour pouvoir identifier une source ponctuelle on dessine l'intersection de plusieurs cones. 
La resolution de ce detecteur pour des photons de 3 MeV est de l'ordre de 1°. 

La tableau 3.4 donne un resume des caracteristiques de trois detecteurs : COS-B 
(mission de l'ESA), EGRET (voir plus haut) et GLAST (mission future de la NASA). La 
mise en service de ce dernier est prevue pour 2005. 

3.4.2 Les detecteurs au sol 

A cause de la loi de puissance du spectre, les flux de gamma au niveau du sol deviennent 
tres petits pour des energies de l'ordre du TeV et au-dela. La nebuleuse du Crabe, par 
exemple, nous donne pour un seuil de 1 TeV un flux de 2 10~ 7 gamma par m 2 et par 
seconde. C'est un flux non mesurable avec un detecteur dans l'espace qui fait 1 m 2 . La 
largeur effective de detection doit etre au minimum de l'ordre de 10 4 m 2 , ce qui conduit a 
mettre en oeuvre des ensembles de compteurs installes au sol. 

Les detecteurs EAS (Extensive Air Shower) et ACT (Atmospheric Cherenkov Telescope) 
font partie de cette categorie d'installations. lis utilisent l'atmosphere comme convertis- 
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COS-B 
1975-1982 


EGRET 
1991-2000 


GLAST 
2005 : 


Domaine E [GeV] 


0.05-5 


0.02-30 


0.01-300 


Detecteur e ± 


chambers a etincelles 


chambers a etincelles 


convert. /Si 


Calorimetre 


Csl 


Nal(Tl) 


Csl, ... 


Resolutions a E=lGeV : 
angulaire 
crE/E 


I 

100% 


1.5 

9% 


0.25 
4.4% 


Champ de vue [steradian] 


0.127T 


0.207T 


1.87T 


Surface de detection 


50 cm 2 


1200 cm 2 


8000 cm 2 


Sensibilite au flux(E>lGeV) 

[cm~ 2 s _1 ] 


10 ~ 6 


10" 8 


10 -io 



Tab. 3.4 - Caracteristiques de detecteurs dans l'espace 





ACT 


EAS 


Seuil en energie 


C(10GeV) 


10 TeV 


Resolution angulaire 


0.1° 


1° 


Champ de vue [steradian] 


10~ 3 (1 source) 


2 


Fraction de temps utile 


~10%(pas de Lune,...) 


-90% 


Rejection hadronique 


100-10000 


tres limitee 



Tab. 3.5 - Caracteristiques de detecteurs au sol 



seur et calorimetre donnant lieu a la formation de gerbes electrophotoniques. 

Les detecteurs EAS mesurent les residus de ces gerbes au sol. De ce fait ils ont un seuil 
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de l'ordre de 10 TeV. Les detecteurs ACT utilisent la lumiere de Cherenkov engendree 
par les particules chargees de ces gerbes et collectee par un ou plusieurs telescopes. Pour 
ce type de detecteurs le seuil est de l'ordre de 10 GeV. 

En general la detection des gamma de grande energie est en competition avec celle de 
hadrons du rayonnement cosmique. A partir de 100 GeV ce bruit de fond hadronique est 
1000 fois plus eleve que le signal des gamma recherche. 

Le tableau 3.5 donne quelques parametres comparatifs de ces detecteurs au sol. 

Les telescopes a lumiere de Cherenkov. 

La lumiere de Cherenkov est emise quand la particule chargee a une vitesse qui depasse 
c/n, ou n est l'indice de refraction du milieu. Pour des electrons de la gerbe e.m. 

£ mi n = m e 7min = 0.511 7min MeV (3.34) 

Dans le cas ou le milieu est l'atmosphere terrestre, on a : 

E miQ « 0.511 x 2(T 1/2 (3.35) 

avec : 

5 = n- 1«2.9 10- 4 e^(-^A_) 
h = hauteur dans l'atmosphere 

Cela donne un seuil de 25 MeV a STP (Standard Temperature and Pressure) et de 
35 MeV a 7.5 km. 

L'angle d'emission du photon de Cherenkov par rapport a la direction de propagation 
de l'electron est au maximum de : 

#max = cos _1 (l/n) (3.36) 

A STP, on obtient max « 1.3°. 

L'intensite d'emission par unite de parcours est : 

A partir de ces relations et d'un bon programme de Monte Carlo qui simule des gerbes 
e.m., on peut determiner la distribution au sol des photons de Cherenkov. Des programmes 
comme CORSIKA permettent de simuler des gerbes de grande energie (y compris des 
gerbes hadroniques) . 

Les resultats sont les suivants : 

1) la lumiere issue d'une gerbe e.m. est assez bien focalisee autour de l'axe de la gerbe ; 

2) le nombre de photons de Cherenkov est proportionnel a l'energie de la gerbe ; 

3) la dispersion temporelle d'arrivee de ces photons est petite (< 5 ns) ; 

4) la surface au sol irradiee par ces photons est C(10 4 ) m 2 . 
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Fig. 3.18 - Simulation de remission Cheren- Fig. 3.19 - Simulation de remission Cheren- 
kov par la gerbe electrophotonique engendree kov par la gerbe hadronique engendree par 
par un gamma de 100 GeV. L'image au sol un proton de 500 GeV. 
est representee dans le diagramme biparame- 
trique. 

La figure 3.18 montre la simulation par CORSIKA d'une gerbe e.m. de 100 GeV a 
incidence verticale. La distribution annulaire des photons Cherenkov est due au fait que 
Tangle d'emission est plus large en fin de gerbe qu'au sommet. Au niveau du sol, les 
trajectoires des photons coincident accidentellement. 

A titre de comparaison, la figure 3.19 est le resultat que l'on obtient pour un proton 
incident de 500 GeV, avec plusieurs sous-gerbes. Cela nous suggere que la mesure de la dis- 
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Fig. 3.20 - Principe de 1'IACT (le gros telescope), et du sampling ACT (l'ensemble des 
capteurs). La gerbe e.m. produit un front d'onde de photons Cherenkov. 




FlG. 3.21 - Le reflecteur op- FlG. 3.22 - Image d'une gerbe produite par un gamma 
tique de 10 m. de TObservatoire de 700 GeV, simulee par CORSIKA. Chaque point repre- 
Whipple. sente un photon. La region limitee par le cercle represente 

quelques degres du ciel. La camera a normalement une 
sensibilite finie et un nombre limite de pixels, de l'ordre 
de 100. 
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tribution laterale de la gerbe fournit un moyen de rejection de la composante hadronique 
primaire. 

II existe deux sortes de ACT : le Imaging ACT (IACT) et le sampling ACT. Le principe 
est represents dans la figure 3.20. 

Dans le cas du sampling ACT, plusieurs detecteurs (au moins 10 stations) observent 
une portion du ciel et enregistrent les distributions temporelle et spatiale des photons. 
Des exemples de ce type sont : ASGAT, THEMISTOCLE et HEGRA. D'anciens fours 
solaires tels que CELESTE ont ete convertis : chaque miroir focalise l'image du ciel sur un 
photomultiplicateur qui se trouve a l'emplacement du four. (Notes : ASGAT, THEMIS- 
TOCLE et CELESTE sont des installations situes sur l'ancienne centrale solaire THEMIS 
dans les Pyrenees Orientales, HEGRA est decrit dans le chapitre 2). 

Dans le cas des IACT, on a essentiellement un telescope comprenant une camera tres 
sensible placee au foyer. Cette camera est constitute en general d'un ensemble (jusqu'a 
une centaine) de photomultiplicateurs. Le prototype en est le reflecteur optique de l'obser- 
vatoire Whipple (figure 3.21). Notons que l'installation HEGRA combine les techniques 
d'imagerie et d'echantillonage. L'image sur le plan focal (figure 3.22) est une petite "come- 
te" qui pointe vers le centre si le telescope est dirige exactement vers la source ponctuelle. 
Sinon, on peut estimer l'erreur a partir d'une parametrisation de l'image (voir Particle 
de Hillas dans Ap. J. 342 (1989) 379-395). Hillas introduit les parametres indiques a la 
figure 3.23. 



Length 




Width 



Fig. 3.23 - Les parametres de Hillas. 

En principe, la detection de deux gerbes permet de determiner la direction de la 
source en tragant l'intersection des axes des deux "cometes". Naturellement la detection 
de plusieurs gerbes est necessaire pour obtenir une bonne resolution angulaire. 

Une amelioration de la technique est l'imagerie stereo, avec au moins deux telescopes. 
Dans ce cas, la superposition des deux images nous donne pour chaque evenement la 
direction de la gerbe (figure 3.24). 

II est tres important d'avoir une bonne resolution angulaire si Ton veut combattre effi- 
cacement le bruit de fond. Ce dernier est du essentiellement au fond cosmique hadronique. 
Si Ton suppose un spectre en loi de puissance E~ 7 (7=1.69 pour la nebuleuse du Crabe et 
1.0 pour MRK421), on peut calculer le rapport Signal/Bruit pour un telescope avec une 
resolution angulaire 89, un seuil d'energie E f , une surface de detection effective au sol A, 
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Fig. 3.24 - Principe de l'imagerie stereo. 



un temps de mesure T et un pouvoir de rejection du bruit cosmique r : 

S_ (E t ) - 85 -'' f~7T 
N a 59 V 1 - r 

Done, si Ton fixe A et T, il faut augmenter r et ameliorer la resolution angulaire 59 
pour augmenter le rapport Dans une certaine mesure la baisse du seuil en energie 
contribue aussi a ameliorer le rapport. La technique stereo devrait permettre d'atteindre 
une resolution d'environ 1/3 de mrad pour une gerbe de 1 TeV. La figure 3.25 montre un 
projet 1 de IACT stereo compose de deux detecteurs a haute resolution (des CCD) et un 
detecteur de trigger central (des PM). 

Le trigger doit etre capable d'extraire un signal estime a 10 photo-electrons (p.e. 2 ) 
par m 2 et par TeV dans un bruit de fond diffus dans la region UV (et les nuits sans lune) 
de 0.45 p.e. ns" 1 m" 2 TeV" 1 . 

Cela est faisable car le signal arrive en moins de 5 ns. De plus il n'occupe qu'une 
portion de l'ordre de 0.1° x 0.3 C tandis que le bruit est isotrope. 

Les figures 3.26 donnent la distribution temporelle de l'arrivee des photons d'une gerbe 
de 100 GeV en fonction de la hauteur a laquelle le photon est genere, ainsi que la dispersion 
temporelle de ces photons. Cette dispersion est effectivement de l'ordre de la ns. Un autre 
point est a relever : la possibility de discriminer gamma et rayons cosmiques hadroniques. 
La topologie de la gerbe nous vient en aide, comme le montre la figure 3.27 qui represente 
le parametre "width" pour des gamma de 100 GeV compares a des protons de 1 TeV (le 
rendement en photons de Cherenkov est environ 10 fois plus faible pour des protons que 
pour des gerbes electrophotoniques). 

^AMT, G. Auriemma et al., Astroparticle Phys. 3, (1995), 215. C. Beyeler, Travail de diplome, IPN, 
Universite Lausanne (1997). 

2 Un p.e. est un electron emis par effet photoelectrique 
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Fig. 3.25 - Stereo IACT a haute resolution, avec un detecteur de trigger au centre (projet 
GAMT). 
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Fig. 3.26 - A gauche, temps d'arrivee des photons Cherenkov au sol en fonction de 
l'altitude du point de creation, pour une gerbe e.m. de 100 GeV. A droite, dispersion 
temporelle. 
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Fig. 3.27 - Le parametre "Width" pour des gamma de 100 GeV et des protons de 1 TeV 
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3.5 Le catalogue de EGRET 



Third RGRFT Catalog 

E>IOOMeV 
+90 




-90 

* Active Galactic Nuclei M Pulsars 

• Unidentified EGRET Sources A Large Mage Ian Cloud (LMC) 

O Solar Flare 



Fig. 3.28 - Le catalogue de EGRET 

Dans le catalogue de 1997, le seuil etait fixe a 100 MeV et une coupure avait ete fixee 
a 5cr sur le bruit de fond. Le resultat est donne dans la figure 3.28. 

On a, entre autre, un ensemble d'environ 50 sources que Ton a associe a des objets 
avec z G [0.03,2.28] (voir eq. 1.8), done probablement extragalactiques. Toutefois a cause 
de la faible resolution angulaire l'association n'est pas aisee. 

3.6 Les sources de gamma multi-GeV. 
3.6.1 Sources galactiques 

La figure 3.29 montre le spectre du Crabe. II s'agit done de remission d'une nebuleuse 
qui est excitee par son pulsar. 

La partie E< 10 10 eV est essentiellement due a la radiation synchrotron, tandis qu'en 
dessus on s'attend a une contribution de l'effet Compton Inverse. Les premieres mesures 
par des ACT confirment les predictions. Dans la figure on a indique aussi la region de 
sensibilite de GLAST et de HESS (un IACT propose par Heidelberg). 
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Fig. 3.29 - Emission de la Nebuleuse du Crabe. 



D'autres signaux sont attendus des SN de type II ou de nebuleuses excitees par des 
sources externes. 

3.6.2 Sources extragalactiques 

Active Galactic Nuclei (AGN) : les "blazars" MRK421 et 450 ont ete detectes a l'ob- 
servatoire Whipple ainsi que par le detecteur HEGRA (voir figure 2.7), avec un seuil de 1 
TeV. Un seuil de 100 GeV ou moins devrait permettre l'observation de plusieurs objets. 

La presence d'evenements explosifs a ete constatee dans MRK421 dans la region du 
keV par ASCA et dans la region du TeV a l'observatoire Whipple (voir figure 3.30). Cela 
souligne la necessite d'avoir une bonne sensibilite sur des intervalles de moins d'une heure. 

Amas de Galaxies : les RC cosmologiques (generes pendant la formation des galaxies) 
devraient produire des pions neutres, done des photons par collision avec les nuages. Ces 
photons devraient etre assez difficilement visibles dans les ACT car la dimension angulaire 
de la source depasse le degre. 

Starbust Galaxies : l'intensite des RC dans le noyau de M82 est estimee etre 100 fois 
plus grande que dans notre galaxie. Cet objet est a 3.2 Mpc de nous, done "ponctuel" et 
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Fig. 3.30 - Evenement explosif detecte a 
l'observatoire Whipple. 
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Fig. 3.31 - La valeur du cutoff attendu sui- 
vant le moment ou les galaxies se sont for- 
mees. Ici l'age est donne par l'indice de "red 
shift" z. 
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Fig. 3.32 - Formation de paires dans le halo d'un objet cosmologique. 



pourrait etre etudie par des ACT. 

3.6.3 Cosmologie observationnelle. 

Absorption photon-photon. Si Ton suppose connu le spectre "original" de photons 
de haute energie, N Q (E), l'interaction avec le fond e.m. sur une distance d en attenue 
Tintensite : 

N(E)=N (E)exp[-r(E,d)] 

Le "cutoff" E* est la valeur de E pour une attenuation 1/e, c. a d. pour r(E*,d) 
= 1. Si Ton peut determiner E* pour plusieurs sources, on peut ensuite comparer les 
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resultats de mesure a des modeles qui contiennent l'age des galaxies comme parametre 
(voir figure 3.31). 

Formation de paires dans le halo d'un objet cosmologique : les gamma de 
la source centrale ne sont pas directement visibles car ils interagissent avec les fonds IR 
(Infra-Rouge) et Optique locaux. Par contre ils produisent des paires e + e~ qui donnent 
lieu a du Compton inverse sur la radiation de fond (voir figure 3.32). On peut montrer 
que cela constitue une "bougie standard" et pourrait amener a un calcul independant de 
la constante de Hubble. 

3.7 Geminga 

C'est une source tres intense qui a ete identifiee a un pulsar avec une periode P=0.237 s 
et un champ magnetique a la surface de 1.6 10 18 T. C'est la seconde source de gamma 
du ciel en importance. On n'a pas pu trouver une contrepartie d'ondes radio, malgre de 
longues recherches. Geminga (en dialecte milanais, signifie "il n'est pas la") fut observee 
par le satellite SAS-2 et etudiee en detail par COS-B. Par la suite, une source X pulsee 
detectee par ROSAT lui a ete associee. La periodicite a aussi ete determinee a l'aide des 
donnees de EGRET, COS-B et SAS-2. 

3.8 Gamma Ray Burst (GRB) 

lis avaient ete observes par les satellites de la serie Vela, qui surveillaient les explosions 
atomiques. L'information fut cachee pendant 13 ans pour des raisons de secrets militaires. 

Le detecteur BATSE a montre que ces evenements explosifs, de courte duree (< 1 sec), 
peuplent de fagon uniforme le ciel. Cela nous fait penser a une origine extragalactique. 

Toutefois, la preuve de leur origine cosmologique ne pouvait etre apportee que par 
l'identification de la source et une mesure du "red shift". En 1997 cela fut realise par 
le satellite X BeppoSAX qui permettait de localiser la position de la source dans une 
region de quelques minutes d'arc. Le 28 02 1997 BeppoSAX signala un evenement qui 
fut reconnu par le telescope Hubble et par le telescope Keck (Sur le Mauna Kea, ile 
d'Hawaii). Ces derniers observerent une nebuleuse de luminosite decroissante. Des lignes 
d'absorption du magnesium ont permis d'estimer un z=0.835, ce qui est un indice de 
l'origine extragalactique de GRB970228. Depuis lors plusieurs objets ont ete correles aux 
observations de BeppoSAX, avec un z maximal de 3.4. 

Quelle est 1'origine des GRB ? 

A la difference des SN, ces evenements sont tres courts, entre une fraction de seconde 
et quelques secondes (figure 3.33). II pourrait s'agir de rejection par une etoile d'une 
couche tres mince de particules relativistes, a la suite d'un choc. La duree tres courte 
que Ton observe dans les diagrammes temporels pourrait etre un effet associe a la vitesse 
ultra-relativiste de l'enveloppe, avec 7 = 0(1000). 

Supposons par exemple qu'un flash 71 ait lieu a t=0 (figure 3.34). A ce moment on 
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Fig. 3.33 - Echantillon de GRB observes par BeppoSAX. 



imagine que la couche soit ejectee avec une vitesse v voisine de c. Alors 

1 1 



7 2 = 



1-/9 2 2(1-/3) 
1 



v = c(3 ph c 1 



et 



2 7 2 

Supposons ensuite qu'un second flash 72 soit genere par la rencontre de la couche avec 
un materiau situe a distance R. Pour un observateur plage au centre, le temps necessaire 
a la lumiere du premier flash pour arriver en R est t=R/c, done la difference temporelle 
entre les deux flashes observee sur terre sera : 



^^terre 



eft c 



R 



t 



2c7 2 2 7 2 



Si Ton considere un temps propre d'un mois, on obtient environ 1 s, pour 7=1000. 

Meme si l'on ne connait pas le mecanisme en detail, on peut considerer qu'on n'est 
pas trop loin de la verite. II a ete suggere que les GRB proviennent de la collision d'etoiles 
a neutrons. Le taux attendu serait, dans ce cas, d' un evenement par 10 6 annees et par 
galaxie, ce qui correspond au taux observe. 
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Temps = t 



Fig. 3.34 - Diagramme pour expliquer la structure temporelle des GRB. 



II a aussi ete suggere que les GRB proviennent de restes de supernovae. Le coeur de 
la supernova (l'etoile a neutrons du centre) emettrait des jets bipolaires (comme c'est le 
cas des noyaux des galaxies actives). Les particules du jet peuvent frapper la couche en 
expansion produite par l'explosion de la SN et la rechauffer localement. Le mouvement 
fortement relativiste transforme le spectre quasi-thermique en une impulsion gamma col- 
limatee. L'hypothese que les GRB et les Supernovae sont des phenomenes lies, semble 
etre confirmee par l'observation de GRB030329 et de l'identification subsequente de la 
Supernova 2003dh (arXiv :astro-ph/0304173, 21-5-2003). 
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4 Neutrini 
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ANNEXES 



A Interaction des particules avec la matiere. 
A.l Photons et electrons. 

On peut decrire phenomenologiquement l'interaction de photons qui traversent une 
epaisseur dx de matiere homogene par l'absoption d'une fraction : 

dN 

— = -pdx 

le signe — indique que les dN sont perdus, \i est le coefficient d'absorption, N le 
nombre de photons a l'entree. 




WfV 1 <[■¥ 1 Hr:V 1 IfliV nfKMIrV 10 c v lkc v IMcV Ll4fV lOIHitV 

Energie du photon Energie du photon 



Fig. A.l - Sections efficaces d'absorption des photons en fonction de leur energie, a) dans 
le carbone ; b) dans le plomb. 



reintegration donne la decroissance exponentielle en fonction de l'epaisseur totale : 

N — N exp(-fix) 

ji peut etre calcule theoriquement a partir de la section efficace a d'absoption, qui repre- 
sente la probability pour un photon d'interagir avec le milieu. Trois processus dominent 
(voir figure A.l) : 

• 1) r : l'effet photoelectrique : 7 + A — > A + + e~, dominant pour E 1 < 500 KeV. 
Les discontinuites dans la courbe cr p . e . correspondent a des transitions entre couches 
electroniques de l'atome ; 

• 2) l'effet Compton (o"c m P ton) : 7 + e _ — > 7 + e~ et la diffusion coherente sur les 
atomes (a Ray i eig h) : 7 + A -> 7 + A ; 

• 3) la production de paires dans le champ nucleaire (k„ mc ) et dans celui des e~(/t e ), 
dominante pour E y > 2 MeV. 
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De plus, on peut observer des reactions photonucleaires avec les noyaux qui sont suivies 
en general par remission de nucleons (oph.n.)- 

La section efficace de l'effet photoelectrique est proportionnelle a Z b pour E 1 > m e c 2 . 
La section efficace de creation de paires est proportionnelle a Z 2 . Le nombre total d'elec- 
trons disponibles pour effectuer une diffusion Compton est de l'ordre de Z/atome. 

En conclusion, l'effet global doit etre fonction du Z (moyen) de la substance. 




Fig. A. 2 - Perte d'energie fractionnaire par unite de longueur de radiation en fonction 
de l'energie des electrons (positrons). 

Un electron qui traverse la matiere perd de l'energie par ionisation comme toute 
particule chargee (figure A. 2). Pour des energies de plus de 10 MeV, le processus dominant 
est celui de Bremsstrahlung : e — > 7 + e qui a lieu dans le champ coulombien d'un noyau. 
L'intervention de ce champ est indispensable pour assurer la conservation de (E,p); il 
s'agit en effet d'une diffusion 7 + e — > 7 + e, ou le premier photon est virtuel et represente 
le quantum d'energie e.m. echangee entre le noyau et l'electron. 

La perte moyenne d'energie par collision depend du Z 2 de la substance. On trouve que 
l'on peut exprimer la perte d'energie par rayonnement sous la forme : 

dE _ E 
dx X 

valable pour E » m e c 2 . X est appelee la longueur de radiation. L'energie moyenne 
des electrons a la profondeur x est donnee approximativement par : 

< E(x) >= E exp(-x/X ) 
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Les fluctuations statistiques sont importantes etant donne la nature stochastique de l'in- 
teraction (theories de Landau et Vavilov). Done, pour x = X l'energie moyenne est 
reduite d'un facteur 1/e. 

La longueur de radiation X depend du materiel traverse ; cette dependance est donnee 
par une expression du type : 



valable a ~20% pour Z > 13. 

Par exemple dans le Fe : 180 • 56/26 2 = 14.9 g/cm 2 =14.9/7.87=1.89 cm, a comparer 
avec 1.76 de la table. 

Si maintenant on revient a la production de photons par rayonnement on trouve que 
la section efflcace est aussi proportionnelle a A/Z 2 . On peut done utiliser la longueur de 
radiation pour caracteriser ce processus. La probabilite de conversion en paire e + e~ est 
donnee par : 



Sur une distance X , la probabilite de conversion d'un photon en un couple e + e~ vaut 
environ 1 — e~ 8 = 0.55, a peu pres constante a partir de 1 GeV. 

En conclusion, a haute energie ou la perte d'energie par radiation et la creation de 
paires dominent, l'interaction des electrons et des photons dans la matiere ne depend que 
de la valeur de X . On peut ainsi caracteriser la topologie des gerbes e.m. initiees par une 
de ces particules de fagon independante de la substance. 

A. 2 Gerbes electromagnetiques. 

Elles sont generees par l'interaction de gamma ou d'electrons energetiques avec la 
matiere (interactions e.m.). Comme precedemment on est amene a normaliser les distances 
en terme d' X la longueur de radiation. Dans ce cas, on obtient pour les electrons : 



De plus il y a une valeur de l'energie de l'electron E c (E critique) en dessous de laquelle 
la perte d'energie par collision (ionisation du milieu) domine la perte par radiation : 



valable a 10% pres pour Z > 13. 

Notons que certains auteurs designent par e l'energie critique et par E c l'energie du 
seuil de detection. 

Un modele simple de generation d'une gerbe e.m. se trouve dans le livre de W. Heitler 3 . 
Un photon incident d'energie E traversant une epaisseur X Q a une prob. = 0.55 de se 
convertir en une paire e + e~. Chaque particule generee a en moyenne une energie E /2 
(figure A. 3). Apres une distance X , l'energie de chaque electron est reduite par un facteur 

3 W. Heitler : The Quantum Theory of Radiation. 3d edition. Oxford Press (1954). 



X w 180A/Z 2 [g/cm 2 ] 




*5 -^radiation — E 



(E > 1 GeV) 



E c = 580/Z MeV 
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1 : 1 Xo 



Fig. A. 3 - Developpement d'une gerbe electromagnetique (modele de W.Heitler). 



er 1 . Une grande partie de cette energie est emise sous forme de photons avec une energie 
comparable a E . On suppose qu'un seul photon est emis avec E = (E /2)/2. On a 
maintenant 4 particules et on est a une profondeur de 2X . Et ainsi de suite. Apres 
tXo, on a 2* particules. Le processus s'arrete quand E /2 t devient plus petit que E c . A 
ce stade, le nombre de particules est maximal : N max = E /E c et la profondeur vaut : 

tmax = hg(E /E c ) log 2. 

Ce modele simple peut etre ameliore. Des programmes de simulation (comme EGS 
"electron gamma shower" de SLAC) sont capables de reproduire les caracteristiques des 
gerbes avec une precision de 1 a 5%. Toutefois des formules simples pour des estimations 
sont disponibles : 

• Position du maximum : t max = \og(E Q / E c ) — t avec t = 1.1 pour des gerbes 
initialisees par des e ± et to = 0.3 pour des gamma. La gerbe e.m. produite par des 
e ± commence plus tot que celle creee par un 7. 

• Le centre de gravite de la gerbe vaut t max + 1.4 et t max +1.7 respectivement pour 
les cas e ± et 7. 

• Un bon detecteur "calorimetrique" doit pouvoir contenir la gerbe e.m. dans sa 
presque totalite. La perte vers l'arriere (mais aussi le "back-scattering" vers l'avant..) 
fluctue autour d'une valeur moyenne et affecte la resolution de la mesure. On trouve 
en moyenne que la gerbe a perdu 98% de son energie dans environ L 98 % = 2.5t max . 
Un bon calorimetre fait 20X au moins. 

• La longueur totale des traces de particules chargees vaut T = E/E c . 

• D'autre part, la dimension laterale de la gerbe est fonction de la profondeur. Physi- 
quement elle est determinee par Tangle moyen de bremsstrahlung p e /m e ) et par 
les effets de diffusion multiple. On trouve que la gerbe est contenue dans un rayon 
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Pm '■ Pm — 21 ^ eV (le rayon de Moliere) jusqu'a la profondeur t max - Integre sur toute 
la profondeur, le rayon necessaire pour collecter 90% de l'energie vaut : Rgo% = ^Pm- 
On a souvent besoin de parametriser le comportement moyen des gerbes. Cela est 
utile pour des routines de "Fast Showering". 

A. 3 Gerbes hadroniques. 

L'interaction e.m. des hadrons charges devient negligeable a haute energie. Les pro- 
cessus forts a partir de quelques GeV/c ont une section efficace totale cr tot de l'ordre de 
20-50 mb, dont une certaine fraction est constitute par de la diffusion elastique (~ 20% 
dans le cas p — p). 

Pour mesurer la probability d'absorption d'un hadron dans la matiere on utilise ge- 
neralement une section efficace qui ne contient pas les contributions elastiques (diffu- 
sion coherente sur les noyaux) ou quasi-elastiques (diffusion sur les nucleons du noyau) : 
°abs — Ctot ~~ °e\ ~~ a q.k\.- P ar exemple on trouve que a a b s ~ 41.2 A - 71 [mb] pour des 
neutrons sur des noyaux de poids atomique A. 

On n'a pas d'approche simple pour decrire une gerbe hadronique. On definit une 
longueur d'absorption : Ao = A/NApcrab& qui permet de "normaliser" la topologie de 
la gerbe (un peu comme X dans la gerbe electromagnetique). Une approximation est par 
exemple : 

A [ cm ] ~ ~ 390 A' 2 / 3 

P 

II existe des programmes de simulation (Geisha, Fluka,...) mais leur precision est reduite. 
Les fluctuations des gerbes hadroniques sont beaucoup plus grandes celles des gerbes e.m., 
ce qui complique la resolution du probleme. 

Une gerbe typique est engendree par la creation de particules rapides par interaction 
forte sur les nucleons de la matiere. Les secondaires sont pour l'essentiel des pions et des 
nucleons. Souvent des n° sont produits, ce qui donne des sous-gerbes e.m. A 10 GeV, 
cela represente 20% de l'energie dissipee : la fraction de n° est donnee a peu pres par 0.1 
log(E [GeV]). Une partie de l'energie est perdue en excitation et fission nucleaire, ce qui 
reduit la part de l'energie "visible". Une autre part de la reduction vient de la production 
de neutrinos non detectes. Dans certaines experiences les muons echappent aussi a la 
detection. 

Rappelons quelques relations utiles : 

• Le maximum de la gerbe (en unites Ao) : t max = 0.90 + 0.36 \og(E[GeV]). 

• La profondeur correspondant a une deposition de 95% de l'energie est donnee par : 
^95% = ^max + 2.5 E 0A3 . En pratique un calorimetre hadronique devrait avoir une 
epaisseur d'au moins 5A . 

• 95% de l'energie de la gerbe est contenue dans un rayon #95% = lAo- 

Le tableau A. 3.1 donne une certain nombre d'informations concernant les materiaux 
usuellement utilises pour les calorimetres. 
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a) 


a) 


a) 


b) 


a) 




<0 


Materiel 


Z 


A 


P 


dE/dx 


A (cm) 


A (cm) 


A 


E c 


e.m. 








(g/cm 2 ) 


(MeV/cm) 


pour N 


pour 7r 


(cm) 


(MeV) 


Aatt(Ao) 


C 


6 


12.0 


w 1.55 


w 2.76 


49.9 


64.9 


w 27.5 


75.9 




Al 


13 


27.0 


2.70 


4.37 


37.2 


45.8 


8.9 


39.3 


2.7 


A liquide 


18 


40.0 


1.40 


2.11 


80.9 


97.7 


14.0 


29.8 




Fe 


26 


55.9 


7.87 


11.6 


17.1 


18.9 


1.76 


20.5 




Cu 


29 


63.5 


8.96 


12.9 


14.8 


17.2 


1.43 


18.7 


3.0 


Sn 


50 


118.7 


7.31 


9.4 


22.8 


24.7 


1.21 


11.4 


3.5 


W 


74 


183.9 


19.3 


22.6 


10.3 


10.5 


0.35 


7.9 


4.1 


Pb 


82 


207.2 


11.35 


12.8 


18.5 


18.4 


0.56 


7.2 


3.9 


U 


92 


238.0 


18.95 


20.7 


12.0 


11.4 


0.32 


6.6 




Nal 






3.67 


4.84 


41.3 




2.59 


12.5 




Scintillateur plastique 


1.032 


2.03 


68.5 




42.9 


87.1 




Lucite 






w 1.2 


2.32 


65.0 




34.5 


80.0 





Tab. A.l - Parametres de materiaux utilises pour les calorimetres. Note : Ao = longueur 
d'absorption nucleaire, X = longueur de radiation, E c =energie critique, a) De la table 
PDG (1978). b) La section efficace inelastique des pions = 31.2 A 744 (mb) a ete 
utilisee. c) Les valeurs correspondent a des gerbes electromagnetiques de 6 GeV. 
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B Annexe : Equation de transport 
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